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ABSTRACT 

Samples of early-type galaxies show a correlation between stellar velocity dispersion and the stellar 
initial mass function (IMF) as inferred from gravity-sensitive absorption lines in the galaxies’ central 
regions. To search for spatial variations in the IMF, we have observed two early-type galaxies with 
Keck/LRIS and measured radial gradients in the strengths of absorption features from 4000-5500 A 
and 8000-10,000 A. We present spatially resolved measurements of the dwarf-sensitive spectral indices 
Nal (8190 A) and Wing-Ford FeH (9915 A), as well as indices for species of H, C 2 , CN, Mg, Ca, TiO, 
and Fe. Our measurements show a metallicity gradient in both objects, and Mg/Fe consistent with a 
shallow gradient in a-enhancement, matching widely observed trends for massive early-type galaxies. 

The Nal index and the CNi index at 4160 A exhibit significantly steeper gradients, with a break at 
r ~ 0.1 r e ff (r ~ 300 pc). Inside this radius Nal strength increases sharply toward the galaxy center, 
consistent with a rapid central rise in [Na/Fe]. In contrast, the ratio of FeH to Fe index strength 
decreases toward the galaxy center. This behavior cannot be reproduced by a steepening IMF inside 
0.1r e ff if the IMF is a single power law. While gradients in the mass function above ~ 0.4 M@ may 
occur, exceptional care is required to disentangle these IMF variations from the extreme variations 
in individual element abundances near the galaxies’ centers. 


1. INTRODUCTION 

Driven primarily by observations within our own 
Galaxy, the assumption of a universal stellar initial 
mass function (IMF) throughout space and time has 
been employed in numerous studies of galaxy evolution. 
The canonical IMFs have the form of a power-law at 
high masses w ith dN/dm oc m~ a and a = 2.3 above 
m > 0.5 M 0 ( IKroupal l2001l) or m > 1M 0 dChabrierl 
I2003T) . These IMFs then flatten and turn-over with 
a decreasing number of stars at masses below 0.5-1 
M 0 . In an extra-galactic context, key observables are 
influenced by the universal-IMF assumption, such as 
the galaxy mass function, mass-metallicity relation, and 
the correlation between total star formation rate and 


galaxy mass (e.g.. iNoeske et al.ll2007HWuvts et al.ll20TTl: 


Leauthaud et al. 

20121 ISmit et al.1120121: Behroozi et al.1 

20131 Zahid et al. 

20141:iSalmon et al. 20151). Spectral en- 


ergy density (SED) fitting is the most prevalent method 
for determining the stellar masses of intermediate- and 
high-redshift gala xies, and requires an assumed form 
of th e IMF (e.g., iConrov et al.l 120091 : iMarchesini et al.1 
I2009T) . Conseque ntly, claims of IMF variations deserve 
intense scrutiny (|Bastian et al.1120101 lKrumhoI3 l2014fl 
since any systemic IMF variation has broad implications 
for the inferred properties of galaxies and galaxy evolu¬ 
tion. 

In the past five years numerous studies have as¬ 
serted that early-type galaxies with the largest ve¬ 
locity dispersions (<j) reveal a bottom-heavy IMF in 
their old stellar populations: an overabundance of 
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stars wi th m < 1 M & relative t o the canonical 

IMFs of IKroupal (120011 1 or IChabrierl (j2003f l. Method¬ 
ologies used to assess the IMF in these galaxies 
have included (1) examinations of stellar absorption 
featu res dominated by either giant or dwarf stars 


2012 

; Conrov & van Dokkum 

|2012bl: Ferreras et al.il 

2013 

: ILa Barbera et al.l 12015 

: ISpiniello et al.1 120141), 


(2) comparisons of mass-to-light ratios from stel¬ 
lar p opulation synthesis (SPS) and stellar dynamics 
(e.g.. iCappellari et al. 120121 l2013bl : IDutton et al.l 120131 : 
iMcDermid et al.l 12014 1. (3) and com parisons between 
SPS and gravitational lensing (e . g.. iTreu et al.1 [20101: 
Barnabe et al.l 120131 : IPosacki et al.l 120151 : ISpiniello et al l 


201 1L 120121 . I2015bl) . Results with the most divergent 
IMFs show that the inferred slope of the IMF above 
0.1 Mq steepens t o a ~ 3 for the most m assive early- 
type galaxies (e.g. JLa Barbera et al.ll2013l ). Other stud¬ 
ies have found a deviation between the average IMF in 
early-type galaxy samples and the canonincal IMF, with¬ 
out verifying a diff erential trend between different early- 
type galaxies (e.g . ■ lAuger et al.ll2010t IDutton et alJl20i2t 
ISmith et al .1120 1 2t ISmith fe Lucevll2013all . Alternatively, 
a few investigations have identified massive early-type 
galaxies with an IMF similar to the Mil ky Way (e.g., 
ISmith fe Lucevll2013bl: ISmith et al.ll2015al) . Using a dif¬ 
ferent approach. IPeacock et al.l ~ 2014fl examined eight 
nearby galaxies for emission from X-ray binaries, and 
found that the fraction of massive stellar remnants across 
their sample was consistent with a uniform IMF slope 
above 8 M 0 . Thus, it is still debated whether the IMF 
varies, how much the functional form changes, and how 
these variations depend on galaxy properties. 

A limitation of the majority of IMF investigations in 
early-type galaxies is the use of a single spatial aper- 






















































































ture per galaxy. For instance, stacked spectra from 
the Sloan Digital Sky Survey (SDSS) have a fixed ra¬ 
dius on the sky, blending data from less than and 
greater than one galaxy effective radius, r R ff, from differ- 
cnt g alaxies ('e.g.JFerreras et alJl2013t lLa Barbera et all 


120131 : iSoiniclio ct ah 20141). In contrast, t he absor ption¬ 
line studies by Ivan Dokkum fe Conrovl (120121) and 
iConrov fc van Dokkuml (l2012bl) focus on the innermost 
regions of nearby galaxies, with a long-slit aperture of 
r e ff/8. For lensing studies, the Einstein radius matches 
a different physical radius in eac h galaxy, typically be¬ 
tween ~ 0 .3r R ff an d ~ lr P f? (e.g. iKoopmans et all 20091 
Smith et al.ll2015aj). iCappellari et al.1 (12012L 12013b ) and 
M cDermid ct ah a2014h use resolved two-dimensional 


stellar kinematics, extending to 1 r e s for many of the 
galaxies in their ATLAS 30 sample. However, the most 
massive galaxies in ATLAS 30 are not covered out to 1 
r e ff. At best, studies with different spatial footprints offer 
leverage for interpreting the role of IMF gradients within 
individual galaxies, although there are numerous compli¬ 
cations from synthesizing heterogeneous and sometimes 
contradictory results. At worst, these studies all fail to 
distinguish between an IMF that varies only from galaxy 
to galaxy and IMF gradients within single galaxies. 

Theoretical motivation for the presence of IMF gra¬ 
dients within early-type galaxies comes from mod¬ 
els of inside-out growth, wherein massive galaxies 
are built first as compact starbursts and then ac- 


Naab et al.1 2001 2014 

: Hopkins et al.l 1201(1 lOser et al.1 

20121: IShankar et al.l 

with observations o: 
galaxies from redshi: 

20131). This model coincides 

size growth in massive red 
Is ~ 2 to the present (e.g., 

Truiil 

o et al.!l2006t van Dokkum et al.l 2010): Patel et all 

2011 

van clcr Wei et al.l 2014 iVulcani et al.l 20141). and 

with the metal-poor stellar halos of nearbv earlv-tvpe 

galaxies fe.g..lCoccato et al.ll2010: Pastorello et alJ 2014t 


iGreene et al.ll2012L 120151) . If the IMF differs between low- 

and high-er galaxies, then the most massive (and highest 
a) early-type galaxies should naturally exhibit IMF gra¬ 
dients, as their outer regions have been assembled from 
smaller systems. 

A second challenge for IMF investigations in early- 
type galaxies is the competing influence of elemental 
abundance ratios, which can drastically reshape stellar 
absorption features and subtly alter the mass-to-light 
ratio of stars. Some SPS models can vary the abu n¬ 
danc e s of individu a l elements (e.g.IGraves fc Schiavonl 
120081 : IConrov fc van Dokkuml I2012al) . Yet the impacts 
of these abundance variations are woefully entangled 
with one another, and with the effects of IMF varia¬ 
tions and other systematics su ch as isochrone offsets in 
temperature- luminosity space (IGr aves fc Schiayonl l2008l : 
IConrov fc van Dokkuml l2012ai r TSpinicllo ct all l2015al l. 
These same models have identified stellar absorption fea¬ 
tures that are especially sensitive to IMF variations, 
and targeted analyses of those features in observed 
galaxy spectra claim suf ficient leverage to detect IMF 
variations robustly (e.g ., Conrov fc van Dokkuml l2012bl : 
lLa Barbera et alJ 120131 : ISpiniello et al.1 120141) . How¬ 
ever, abundance ratio s are known to va r y within indi¬ 
vidua l galaxies (e.g., Strom et al.l 119761: iTamura et al.1 


20001: iWeiimans et al. 1 20091 iKuntschner et al.1 120101 : 


Greene et a.1.112012L120131120151) . and the few investiga 


tions of IMF gradients in individual systems have of¬ 
fered scant analysis of single-elem ent abundance varia¬ 
tions (IMartfn-Navarro ct al .1 120 15 al fbUcil . 

Herein we examine two early-type galaxies, NGC 1023 
and NGC 2974, for IMF and abundance gradients. For 
each object, we use long-slit spectra to probe spatial 
scales from ~ 100 pc to a few kpc (0.03 r e g to ~ 1 r e g). 
We present gradients in a selected set of stellar absorp¬ 
tion line indices and qualitatively interpret their connec¬ 
tion to stellar population gradients. We have paid spe¬ 
cial attention to the sodium doublet at 8190 A (here¬ 
after Nal) and the Wing-Ford iron hydride feature at 
9915 A (hereafter FeH), both of which are sensitive to 
the number density of cool dwarf stars. Our spectra also 
cover the giant-sensitive calcium triplet at 8500-8660 A 
(hereafter Call), temperature-sensitive TiO features, nu¬ 
merous features of individual atomic species, and several 
Balmer lines._ 

Recently, IMartfn-Navarro et al.l (l2015al TbT) analyzed 
long-slit spectra of four nearby early-type galaxies and 
reported IMF gradients in two ellipticals with large cen¬ 
tral <7 (270-300 kms -1 ) and a mild gradient in the 
compact, high-cr galaxy NGC 1277. They inferred a 
uniform IMF slope in an elliptical galaxy with a ss 
100kms -1 . The two objects presented herein have cen¬ 
tral a ~ 210-250 kms -1 , and our set of IMF-sensitive 
spectral indices has lit t le over lap with those analyzed by 
IMartfn-Navarro et all (|2015al fbh . Although our spectral 
coverage excludes the TiOi and TiC >2 features at 5960 A 
and 6230 A, we perform much more rigorous analysis of 
Nal and are among the first to present spatially resolved 
measurements of FeH. 

We summarize the basic properties of NGC 1023 and 
NGC 2974 in Table [T| NGC 1023 is an SB0 galaxy, and 
NGC 2974 is an E4, as classified by the NASA/IPAC 
Extragalactic Databas e. Both galaxies are fast rotators 
(lEmsellem et al.l 120 111). We selected b oth objects from 
the sample of Ivan Dokkum fc Conrovl (|2012D . who ob¬ 
served the central r e ff/8 of 34 early-type galaxies with 
identical spectral coverage to our investigation. Based 
on full- spectrum fitting to st ellar population models, 
IConrov fc van Dokkuml (l2012bl) measured a stellar mass- 
to-light ratio, Tk, of 1.53Y;yMW in NGC 1023, and 
1.45 Tx,mw in NGC 2974, where T*-,mw is the stellar 
mass-to-light ratio of the Milky Way disk. They infer 
t hat the central IMF in both galaxies is consistent with 
a lSalpeteH (119551 1 power-law (a = 2.35) extending down 
to 0.1 Mq and significantly more bottom-heavy than the 
Kroupa or Chabrier forms. 

This paper is organized as follows. We summarize our 
long-slit observations in (}2]and our data reduction meth¬ 
ods in U3 In [2] we describe how absorption line indices 
are determined from our spectra, including random and 
systematic errors. We present radial trends for 13 se¬ 
lected line indices in tj5j In iJ7]we offer a qualitative inter¬ 
pretation of the underlying stellar population trends in 
NGC 1023 and NGC 2974, with rigorous stellar popula¬ 
tion modeling deferred for future work. Our conclusions 
are briefly summarized in fj8] Our two appendices con¬ 
tain detailed examinations of systematic errors in line 
index measurements, with particular focus on the Nal 
feature. 
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Table 1 

Targets and Observations 


Galaxy 

D 

(Mpc) 

M k 

r ef£ 

(") 

cr c 

(kms -1 ) 

[Fe/H] c 
dex 

[Mg/Fe] c 

dex 

[Z/H] c 

dex 

[a/Fe] c 

dex 

Age c 

Gyr 

Integration time 

to 

Slit PA 

(°) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

NGC 1023 

11.1 

-24.01 

47.8 

217 

-0.01 

+0.18 

+0.09 

+0.19 

11.7 

13 x 600 

85 

NGC 2974 

20.9 

-23.62 

38.0 

247 

-0.06 

+0.20 

+0.11 

+0.29 

8.9 

14 x 600 

45 


Notes: We adopt dist ances, A-band absolute magnitudes, and effective radii (columns 2-4) fro m the ATLAS 30 survey 
i Cappellari et, all 120111 1. Central velocity dispersions, [Fe/H], and [Mg/Fe] (columns 5-7) are from IConrov fc va n Dokkum 
i 2012bj}_and appro ximate a circular aperture with radius r e ff/8. Central [Z/H], [a/Fe], and stellar ages (columns 8-10) are from 
iKuntschner et all (1201(1 1. summing SAURON integral-field data ove r a circular aperture of r e s/8. Using the new data herein 
and stellar population models from IConrov fe van Dokkuml (l2012al i. we estimate that both galaxies are > 10 Gyr old at all 
radii. 


2. OBSERVATIONS 

We observed NGC 1023 and NGC 2974 in Decem¬ 
ber 2013 using LRIS on t he Keck 1 telescope (|Oke et al.l 
119951 iRockosi et aI1l2010D . in long-slit mode. LRIS com¬ 
prises a red arm and a blue arm. For the red side, we 
used the 600 line mm -1 grating spanning 7500-10,800 A, 
sampled at ~ 0.79 A per pixel. On the blue side we 
used the 600 line mm -1 grism spanning 3100-5560 A, 
sampled at ss 0.62 A per pixel. For this instrumental 
setup we measured a spectral resolution of « 3.0 A full 
width at half-maximum (FWHM) for both the blue arm 
and the red arm. Our slit was 0.7" wide and spanned 
a length of 175". The LRIS field of view is sampled by 
two detectors on each spectrograph arm, resulting in a 
coverage gap near the center of the slit. For the data pre¬ 
sented herein, the gap spans 35" on the red side and 14" 
on the blue side. On the LRIS detectors we employed 
a pixel scale of 0.135" in the spatial direction. Seeing 
on the night of observations was between 1.0" and 1.5" 
FWHM. 

We placed our slit along the major axis of each galaxy 
with the galaxy center slightly offset from the chip gap. 
For a single exposure our spatial coverage extended from 
« 60" on one side of the galaxy to « 115" on the op¬ 
posite side. Images were dithered between two positions 
that straddled the gap to achieve symmetric coverage of 
each galaxy. Additionally, a sky field 4' away was ob¬ 
served after every two to four science exposures. Total 
integration times for each target are included in Tabic [TJ 
We completed more than two hours of integration time 
on each galaxy in order to obtain signal-to-noise ratios 
( S/N ) ~ 100 per pixel in spatial bins near lr e g. The 
final spectra for NGC 1023 meet this criterion out to 
1.6r e ff, for all features except the faint FeH band near 
1 /im. For NGC 2974 the outermost bin with sufficient 
S/N spans 0.5-0.8 r e g. 

Although our instrument setup matched that of 
Ivan Dokkum fe Conrovl (I2012D . our observing pattern 
was selected to maximize spatial coverage and measure 
gradients i n absorption lin e depth s. In contrast, obser¬ 
vations bv Ivan Dokkum fc Conrovl (12012 1 aligned the slit 
along the minor axis of each galaxy and used the far edges 
for in-frame sky subtraction. As a result, their previous 
analysis was restricted to the central few arc seconds of 
each galaxy, corresponding to an aperture of r e s/8. 


3. DATA PROCESSING AND ANALYSIS 

Our data reduction procedures largely follow the 
template described in Ivan Dokkum fe Conrovl (I2012D . 
T he notab l e exception is sk y subtraction. Whereas 
Ivan Dokkum fc Conrovl ( 2012 1 aligned the slit with the 
minor axis of each galaxy and used the far edges of the slit 
for sky subtraction, we wish to extract spectral features 
over the entire slit. To this end we recorded separate, 
non-concurrent sky exposures. On the red side, our sky 
subtraction requires careful wavelength calibration for 
each science frame, and a scaled subtraction procedure 
that adjusts the relative strengths of telluric emission 
line families, which vary on timescales of a few minutes. 
Details of these various calibration steps are described 
below, and representative cleaned spectra are illustrated 
in Figures [Hand [2] 

3.1. Instrumental Calibrations 

We use the lpipe software packag(0 (Daniel Perley, 
private communication) for bias subtraction, flatfield cor¬ 
rection, and cosmic ray cleaning. This package first de¬ 
termines and subtracts bias levels from the overscan re¬ 
gion of each frame. Next it performs flatfielding using 
halogen lamp exposures taken at the end of the observ¬ 
ing night, with the slit and disperser in place. The lpipe 
package coadds individual flatfield exposures, computes 
a boxcar-smoothed response spectrum for each row of the 
coadded flatfield, and divides each row by its matching 
response spectrum to produce a pixel-to-pixel flat. This 
flatfielding correction retains response variations along 
the spatial axis of the CCD, which are calibrated dur¬ 
ing the sky subtraction step ( 33.21) . After flat-fielding we 
split the raw frames into the separate CCD chips for each 
arm and process each chip independently. 

On the red side, we use telluric OH emission lines to 
define a wavelength solution for each frame, using the 
IRAF routines identify, fitcoords, and transform. 
We divide each frame into blocks of 40 rows (5.4") and 
extract a one-dimensional spectrum for each block. We 
fit a fifth-order polynomial to the peaks of the OH lines, 
converting pixels to A. To calibrate for wavelength vari¬ 
ation across each chip, we then fit a two-dimensional 
polynomial solution to the set of peak locations from 

4 The lpipe IDL routines are available at 
http://www.astro.caltech.edu/~dperley/programs/lpipe.html 
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Figure 1. LRIS blue arm spectra extracted from the central and outer regions of NGC 1023 (top, green) and NGC 2974 (bottom, blue), 
and shifted to rest-fra me w avelengths. Each spectrum has been sky-subtracted and flux-calibrated. Emission lines have been removed, as 
described in Appendix I A. 2| The spectra are normalized to the same median value, and constant offsets have been added for clarity. Stellar 
absorption features analyzed herein are labeled. 
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Figure 2. LRIS red arm spectra extracted from the central and outer regions of NGC 1023 (top, green) and NGC 2974 (bottom, blue), 
and shifted to rest-frame wavelengths. Each spectrum has been background-subtracted with a scaled sky spectrum, corrected for telluric 
absorption, and flux-calibrated. The spectra are normalized to the same median value, and constant offsets have been added for clarity. 
Stellar absorption features analyzed herein are labeled. The increased noise at 9200-9700 A corresponds to strong telluric absorption bands. 
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all blocks. The 40-row extraction window permits reli¬ 
able sky line fitting except for one or two blocks near the 
center of each galaxy; these blocks are masked from the 
two-dimensional fit. 

The blue side includes a faint telluric line from NI at 
5199 A but is otherwise devoid of telluric emission fea¬ 
tures. We therefore derive a two-dimensional wavelength 
solution from daytime arc lamp frames, and assume that 
temporal variation in the wavelength solution can be de¬ 
scribed by a constant offset term at all wavelengths, along 
the entire slit. In some frames, the NI line is too faint to 
perform a useful fit. We therefore use the Fe5270 absorp¬ 
tion feature at the galaxy center to measure the relative 
wavelength offsets across a sequence of frames. We an¬ 
chor the frame-to-frame offsets to a single exposure for 
each galaxy, where the NI 5199 feature is strong enough 
to establish the absolute wavelength scale. In order to 
avoid substantial velocity structure when measuring the 
wavelength shifts in Fe5270, we only extract a 1.5" region 
at the galaxy center. We account for each galaxy’s rota¬ 
tion curve once we have extracted final one-dimensional 
spectra for all spatial bins, and before we measure line 
indices (see |J4|). 

3.2. Sky Subtraction and Telluric Correction 

The background spectrum at red wavelengths is domi¬ 
nated by telluric emission lines, whose relative strengths 
vary on timescales comparable to our exposure times. 
Our goal of extracting galaxy light over the full length of 
the LRIS slit prohibits in-frame sky subtraction, and our 
sky frames must be corrected for variation in the line 
strengths. We extract a high-S/N sky spectrum from 
each science frame by applying a large aperture, offset 
from the galaxy center by at least 30". We collapse the 
entire sky frame into a one-dimensiona l spectrum and 
use the Skye or r routine dNoll et al.l I2014T) El to perform 
scaled sky subtraction on the spectrum extracted from 
our science frame. Skycorr adjusts the relative ampli¬ 
tudes of OH and O 2 emission groups in the input sky 
spectrum to best match the input science-frame spec¬ 
trum, and subtracts the rescaled sky spectrum to output 
a clean science spectrum. 

Subtracting the Skycorr output spectrum from the 
initial science-frame spectrum yields a one-dimensional 
“master” sky spectrum for that particular frame. Next 
we use the halogen flats to compute the average response 
function of the CCD chip in the spatial dimension. We 
expand the master sky spectrum as a two-dimensional ar¬ 
ray, scaled by the spatial response function, and subtract 
this array from the science frame. Figure [3] illustrates 
the improvement in sky subtraction after performing this 
scaled sky procedure, relative to direct sky subtraction. 

Residuals from telluric emission lines are our dominant 
source of noise near the FeH, MgI0.88, and TiOO.89 fea¬ 
tures (see Table [2] for definitions). These residuals are 
likely a combination of high shot noise from the bright 
lines, and wavelength calibration errors with magnitude 
< 0.1 A. We have experimented with multiple variants 
of our sky subtraction procedure, including zoomed-in 
wavelength calibration over a narrow wavelength inter- 

Skycorr and SkyCalc are available 

from the European Southern Observatory at 
http: //www.eso.org/sci/software/pipelines/skytools 
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Figure 3. Sky-subtracted spectra extracted from +30 ,/ to +50 // 
with respect to the center of NGC 1023. Top: region near the 
Ca II feature. Bottom: region near the FeH feature, with dark and 
light shaded areas corresponding to the line center and pseudo¬ 
continuum regions. In both panels, the red (upper) spectrum was 
obtained by directly subtracting sky frames from science frames, 
whereas the blue (lower) spectrum is our final spectrum after scaled 
sky subtraction. Even after scaled sky subtraction, pixels cor¬ 
responding to bright sky lines exhibit excess noise. The dashed 
regions in each spectrum indicate areas that we masked before 
computing line indices. 
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val near the FeH band, fitting the wavelength solution 
with different polynomial orders, and applying Skycorr 
over narrower chunks of wavelength space. None of these 
attempts yielded clear improvements. 

In addition to their native pixels, residuals from 
bright sky lines may contaminate adjacent pixels dur¬ 
ing smoothing, which we employ to bring spectra to a 
common a O- To mitigate this, we perform a ver¬ 
sion _of the i terative masking procedure described by 
Ivan Dokkum fe Conrovl (I2012T) . We first construct a 
mean sky spectrum for our observing night and flag the 
pixels corresponding to the brightest sky lines. We then 
mask these pixels at the corresponding wavelengths in 
each galaxy spectrum, such that their values are inter¬ 
polated from nearby good pixels before smoothing. As 
the smoothing introduces undesired correlations between 
good pixels and masked pixels, we return to the un- 
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smoothed spectrum and replace only the masked pixels 
with the smoothed output. We then perform the smooth¬ 
ing again, and iterate the substituting and smoothing 
steps five times. Although pixels near the edges of our 
masked regions are still over-weighted in our final spec¬ 
trum, the iterations serve to distribute the excess weights 
more broadly. 

In Figure [3)o it is evident that the sky masking removes 
a substantial fraction of the FeH line region at 9902- 
9932 A. In Appendix [B] we examine additional variants 
of the FeH index that extend the line region to 9962 
A. These variants contain a larger fraction of unmasked 
pixels and yield similar radial trends in the FeH index 
for r < 0.5r e ff. As a second test, we have smoothed our 
spectra and measured line indices without any masking. 
The resulting index values do not differ significantly from 
the measurements we present below. Only the masked 
data exhibit a marginal upturn in MgI0.88 at large radii 
in NGC 1023 (Figure [5]). Otherwise, the radial trends 
in MgI0.88 and FeH are qualitatively similar for masked 
and unmasked measurements. 

After subtracting the sky emission spectrum, we cor¬ 
rect each red arm frame for telluric absorption. We start 
with th ree transmission spectra from the SkyCalc sky 
model (|Noll et all 120121 1 Jones et~aI112013H 0 a baseline 
spectrum ( T B ) for low airmass and low water vapor, a 
high-airmass spectrum (Ta), and a high water-vapor 
spectrum ( Tw )• After convolving each model spectrum 
to the instrumental resolution of the LRIS red arm, we 
form a grid of linear combinations (Tc): 

Tc =T b + v(Ta ~ T b ) + uj(T w - T B ) + C • (1) 

For each exposure, we extract a spectrum from the cen¬ 
tral region of the galaxy, flatten it over the range 9250- 
9650 A by dividing out a fourth-order polynomial, and 
determine the values of ( 77 , ui, () that best reproduce the 
observed telluric absorption over the wavelength range 
9310-9370 A. We then divide the frames on both CCD 
chips by the corresponding model Tc . At small radii, 
noise from our transmission correction dominates the 
galaxy spectra over the range ~ 9200-9700 A, where 
there are no stellar absorption features of interest. More 
importantly, a band of telluric H 2 O overlaps the red- 
shifted Nal feature in both galaxies, with a particularly 
strong trans missi on dip at 8230 A observed wavelength. 
In Appendix lA.3l we assess a plausible error range for the 
transmission spectra used to correct the Nal feature, and 
the corresponding error e te i in our measurements of the 
Nal line index. We find that variations of 10% in our 
adopted transmission spectrum lead to errors as high as 
13% and 4% in the Nal index for NGC 1023 and NGC 
2974, respectively. These errors are incorporated along 
with other systematic terms in the line index measure¬ 
ments presented below (e.g., in Figures [5] and ©• We 
note that our empirically derive d errors are mu ch larger 
than t he estimate of 0.1-0.2% bv lvan Dokkum fo Conro\l 

mm- 

For our wavelength coverage on the blue side, the sky 
background is dominated by a continuum spectrum, and 
the background level was very stable during our observ¬ 
ing night in December 2013. For the corresponding blue- 
arm data we directly subtract a calibrated sky frame 


from each calibrated science frame. Telluric absorption 
is negligible on the blue side. 

3.3. Position Registration 

For each science exposure, we trace the position of the 
galaxy center by fitting a two-component Gaussian pro¬ 
file to the central « 5" of each galaxy, in each of 10 
wavelength blocks. The resulting trace is interpolated 
to all wavelengths, and one-dimensional spectra are ex¬ 
tracted with sub-pixel spatial precision at each wave¬ 
length, approximating the flux in each pixel to be evenly 
distributed in the spatial dimension. Each science expo¬ 
sure also includes one chip per arm that is offset from 
the galaxy center. We measure the spatial gap between 
the on-center and off-center chip by comparing the RA 
and DEC header keywords for alternating dither posi¬ 
tions to the location of the galaxy center on the corre¬ 
sponding chips. For NGC 1023 and NGC 2974 the gap 
spans 35" on the red arm and 14" on the blue arm. This 
includes a small buffer region where the data frames are 
trimmed during initial calibration steps. To calibrate for 
wavelength-dependence of the trace on each off-center 
frame, we apply the trace from an adjacent exposure 
where the galaxy center is positioned on the correspond¬ 
ing chip. 

3.4. Spectral Response Calibration 

On t he red side, we fo ll ow the flux calibration proce¬ 
dure of Ivan Dokkum fc Conrovl (120121) . which uses halo¬ 
gen flats to measure variations over small wavelength 
scales and calibrates large-scale variation with a white 
dwarf spectrum. We extract a one-dimensional halogen 
spectrum for each chip, apply 10-pixel boxcar smooth¬ 
ing, and divide the result into the extracted spectrum of 
a white dwarf calibration star. After this initial correc¬ 
tion step, the white dwarf spectrum should approximate 
a A -4 power law, with deviations on ~ 100 A scales 
arising from the intrinsic shape of the halogen spectrum. 
We divide the corrected white dwarf spectrum by a A~ 4 
function and fit the residual spectrum with a 10-point 
cubic spline. We then multiply the smoothed halogen 
spectrum for each chip by this smoothed residual pro¬ 
file, flattening the intrinsic source and yielding our final 
response curve. 

On the blue side, the steep intrinsic spectrum of the 
halogen light source renders even preliminary calibration 
ineffective, and so we divide the observed white dwarf 
spectrum on each chip directly by A~ 4 . To remove stel¬ 
lar features and pixel-to-pixel noise, we mask the pixels 
corresponding to Balmer absorption features, interpolate 
from the remaining wavelengths with a cubic spline, and 
perform 100-pixel boxcar smoothing to obtain a final re¬ 
sponse curve. Visual inspection of our blue-side halogen 
spectra reassures us that there is little response variation 
on scales < 50 A. 


4. MEASURING LINE INDICES 

While some investigations of stellar pop ul ations 
empl oy full spectral fitting l e.g.. iCid Fernandes et al.1 


200l iKoleva et al.l 20091 

Kuntschner et al.l 

2010 ; 

Conrov & van Dokkum 2012bl: 

Poclorvanvuk et al.l 

2013: 

Conrov et al. 2014: McDermid et al.ll2015t Posacki et al. 


20151 iWilkinson et al.l l2015f> . line indices or equivalent 
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widths of specific absorption features are useful for 
qualitative interpretation and ca n be readily applied to 
a number of SPS models (e.g.. iTrager et all l2000aUbl: 
Thomas et al.l 120051: iSchiavonl 12007 : Grave s fc Schiavonl 


20081 lLa Barbera et all 2013t iSpiniello et al.l I2014T ). In 


Tabic [5] we list the definitions of 13 line indices discussed 
herein. On the blue side we track promi nent i ndices 
from the Lick/IDS system introduced by iFaber et~ahl 


Tracer et all 

1991. 

Spiniello et al. 

(1201411 


as well as the bTiO index from 
Our wavelength coverage with 
the LRIS blue arm cuts off in the middle of their aTiO 
feature. At near-infrared wavelengths, we adopt line 
i ndex definitions for Na l, Cal l, FeH, TiO, and Mg from 
iConrov fe van Dokkunj (12012ai hereafter CvD12). 

lLa Barbera et ahP 20131) introduced an alternative def¬ 
inition for the Nal index, whose pseudo-continuum re¬ 
gions traced smaller residuals between stellar population 
models and galaxy spectra from SDSS. We describe this 
feature, NalsDSS, in Table[2j In Appendix iBlwe examine 
the radial behavior of NalsDSS and three new variants 
of the FeH index, to test for biases that might arise from 
overlapping absorption features or contamination by sky 
lines. In brief, we find that the choice of Nal or FeH 
definition does not change the essential radial trends we 
observe (SJ5J), or our interpretation (){7|). 

Our spatial binning scheme for each galaxy is designed 
to maximize S/N in each aperture, while retaining the 
ability to compare data near 0.5r e ff and 1.0r e ff. Be¬ 
tween ~ 0.2 r e ff and ~ 1 r e g our bin sizes range from 
5" to 20". At smaller radii, increased surface bright¬ 
ness permits us to employ much finer binning, down to 
seeing-limited scales of 1.0". In each two-dimensional 
science frame we determine an aperture according to the 
spatial registration described in 43.31 and at each wave¬ 
length we measure the mean flux over the corresponding 
rows. Our spatial mean includes 3 -ct clipping at radii 
where the surface brightness profile is sufficiently shal¬ 
low, typically r > 5". Once spectra have been extracted 
from individual science frames, they are coadded via di¬ 
rect summation. 

To assess spatial gradients in line indices, we must com¬ 
pare spectra with the same velocity dispersion a, in a 
common rest frame. To this end, we measure a and the 
radial velocity v for each spatially binned spectrum (4j5l). 
and apply a Gaussian smoothing kernel so as to arti¬ 
ficially raise a to a common value. The kernel width 
is chosen such that each binned spectrum has a final 
dispersion of 230kms _1 for NGC 1023 and 245kms _1 
for NGC 2974, matching our highest measurement of <r 
in each galaxy. We also shift the wavelength grid for 
each spectrum to rest-frame wavelengths, based on our 
measurement of v. Figure [4] illustrates spectra near the 
Nal, Call, and FeH features, after performing kinematic 
calibrations. Our measurements of v and a employ the 
pPXF proc edur e by iCannellari fe Emselleml (120041) . In 
Appendix IA.1I we discuss possible systematic errors in 
measuring kinematics, and the resulting impact on our 
measurements of line indices. 

The blue-arm spectra of NGC 2974 include strong 
emission lines, which must be removed before measur¬ 
ing the depths of nearby absorption features. We fit and 
subtract an emission line component for each blue-arm 
spectrum by including Gaussian emission line profiles for 


H, NI, Oil, and OIII in our list of kinematic templates 
for pPXF. We explore uncertainties and alternative meth¬ 
ods for the emission line fitting in Appendix I A. 2 1 NGC 
1023 shows much subtl er tr aces of emission, which we 
also discus in Appendix IA.2I 

Once our spectra are cleaned and calibrated for kine- 
matic s, we com p ute equivalent widths using the formu¬ 
lae of lWorthev et al.l (11994 Equations 1-3), whereby the 
continuum level is modeled as a straight line connect¬ 
ing the midpoints of the red and blue pseudo-continuum 
bands in Table [2] Our only exception is the TiOO.89 in¬ 
dex, which is expressed as the r atio of b lue to red pseudo¬ 
continuum levels as defined by lCvD12t 

4.1. Random and Systematic Errors 

In principle, noise in a galaxy spectrum can be prop¬ 
agated analytically to compute formal statistical errors 
in an ensuing equivalent width measurement. Yet noise 
from sky subtraction and telluric absorption is not ran¬ 
dom and uncorrelated in our spectra. Instead of formal 
error propagation, we perform a simpler and more empir¬ 
ical estimate of line index measurement errors, as follows. 
We split the science frames for each galaxy into two sub¬ 
sets, with frames from each dither position evenly dis¬ 
tributed. We extract one-dimensional spectra from each 
frame, matching the spatial apertures in [J51 and coadd 
all the spectra within a subset. After pairing data from 
opposite sides of the galaxy and measuring indices as 
described above, we have four index measurements for 
each interval in r. We adopt the mean and standard 
deviation of these measurements as our final equivalent 
width and ltr random error. We then add the random 
errors in quadrature with the total systematic error e sys 
from kinematic uncertainties, emission line removal, and 
telluric absorption. We derive e sys in Appendix [A] and 
summarize its value for each index in Table [3] 

For most of the line indices we have investigated, e sys 
is comparable to the random deviations between our four 
subsets of data, particularly at radii 10" where our 
spectra have very high S/N. Exceptions where random 
variance dominates systematic errors at all radii are the 
bTiO index in NGC 1023, the Call index in NGC 2974, 
and the TiOO.89 and FeH indices in both galaxies. For 
the FeH feature in particular, the noise including sky 
lines exceeds the 3.5% total systematic error by at least 
a factor of two. Our error bars in Figures [5] and [6] below 
include both random and systematic terms. 

The largest systematic effects occur for Nal and are 
discussed extensively in Appendix [A] In particular, we 
note that the sizable error bars for Nal in NGC 1023 
(Figure [5ji) are largely driven by the ~ 10% etei term for 
possible errors in telluric correction. Biases introduced 
by imperfect telluric correction vary gradually with ra¬ 
dius, as rotation in each galaxy’s rotation curve shifts 
the Na l fea ture across the overlapping telluric band (Ap- 
pendix lA.3l) . While this term is important for considering 
the absolute strength of the Nal index, the hypothetical 
telluric bias may be approximated as a uniform offset 
near the center of each galaxy. Therefore, radial trends 
in Nal are more significant than suggested by Figure [5] 
In NGC 2974, e te i is < 4% and is outweighed by other 
terms in the total error budget. When e te i is excluded, 
the remaining systematic error in Nal ranges from 7- 
12 %. 
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Table 2 

Line Index Definitions 


Index 

Ref. 

Line 

(A) 

Blue Pseudo 

(A) 

Red Pseudo 

(A) 

Units 

Dependence 

CNi 

1 

4142.1 - 4177.1 

4080.1 - 4117.6 

4244.1 - 4284.1 

mag 

ct, Nt, 01 , aget 

Ca4227 (CaI0.42) 

1 (2) 

4222.2 - 4234.8 

4211.0 - 4219.8 

4241.0 - 4251.0 

A 

Caf, Ct, aget 

C 2 4668 (C 2 O. 47 ) 

1 (2) 

4634.0 - 4720.2 

4611.5 - 4630.2 

4742.8 - 4756.5 

A 

Ct, aget, Ot 

bTiO 

3 

4758.5 - 4800.0 

4742.8 - 4756.5 

4827.9 - 4847.9 

mag 

T efr t, g t, Mgt, ot, Tit, Ct 

H/3 

1 

4847.9 - 4876.6 

4827.9 - 4847.9 

4876.6 - 4891.6 

A 

aget, Ct, a t, Fet 

Mg b (MgI0.52a) 

1 (2) 

5160.1 - 5192.6 

5142.6 - 5161.4 

5191.4 - 5206.4 

A 

Mgt, aget, Ct, T eS t 

Fe5270 (FeI0.52) 

1 (2) 

5245.6 - 5285.6 

5233.2 - 5248.2 

5285.6 - 5318.2 

A 

Fet, aget 

Fe5335 (FeI0.53) 

1 (2) 

5312.1 - 5352.1 

5304.6 - 5315.9 

5353.4 - 5363.4 

A 

Fet, aget 

Nal (NaI0.82) 

2 

8174.8 - 8202.7 

8167.8 - 8174.8 

8202.7 - 8212.7 

A 

Nat, 9 t, T eS t, ol t, aget 

Nalgngg 

4 

8177.8 - 8197.7 

8140.8 - 8150.8 

8230.7 - 8241.7 

A 

Nat, g t, T eS t, a t, aget 

CaII=Cal+Ca2 +Ca3 







Cal 

2,4,5 

8481.7 - 8510.7 

8471.7 - 8481.7 

8560.6 - 8574.6 

A 

Cat, g t, Nat, T eS t, Mgt, Fet 

Ca2 

2,4,5 

8519.7 - 8559.6 

8471.7 - 8481.7 

8560.6 - 8574.6 

A 

Cat, g t, Nat, T eS t, Mgt, Fet 

Ca3 

2,4,5 

8639.6 - 8679.6 

8616.6 - 8639.6 

8697.6 - 8722.6 

A 

Cat, g t, Nat, T eS t, Mgt, Fet 

MgI0.88 

2 

8799.5 - 8814.5 

8775.0 - 8787.0 

8845.0 - 8855.0 

A 

Mgt, aget, TeS 

TiOO.89 

2 


8832.6 - 8852.6 

8867.6 - 8887.6 

ratio 

T eS t, Ot, Tit, Ct 

FeH (FeH0.99) 

2 

9902.3 - 9932.3 

9852.3 - 9877.3 

9937.3 - 9967.3 

A 

TeS t, g t, Fet 


Notes: The index definitions above use air wavelengths, as do all figures in this paper. The TiOO.89 index is defined as the 
flux ratio for the blue pseudo-continuum divided by the red pseudo-continuum. The last column lists some of the main atomic 
species and stellar atmosphere properties that influence the depth of each feature. For integrated light, T e s reflects the relative 
contribution of warm vs. cool stars, surface gravity g reflects the relative contribution of giants vs. dwarfs, and a reflects the 
combined abundances of a-process elements. Arrows indicate whether the absorption index increases or decreases with respect 
to an increase in the respective quantity. For instance, g ? means the index becomes stronger with increasin g surface gravity ; 
i.e. the feature is dwarf-sensitive. These dependencies have been compiled from several previous studies fe.g. lServen et al.H2005l : 
ISchiavonll2007l . CvD12). Column 2 references are: (1) Trager et al. 1998; (2) CvD12 (3) Spiniello et al. 2014; (4) La Barbera 
et al. 2013; (5) Cenarro et al. 2001. 


Moreover, Ivan Dokkum &; Conrovl (120121) have sug¬ 
gested that the scatter in the relation between Nal 
strength and the center of the Nal+TiO blend provides 
a heuristic upper limit for the total error in the mea¬ 
sured Nal index. In essence, if the only varying quantity 
is the Nal component of the Nal+TiO blend, this will 
drive a tight anti-correlation between the blend center 
and the Nal index strength. Additional scatter in the 
relation reflects a combination of measurement errors in 
the Nal index and blend center, and independent varia¬ 
tion of the TiO component. We have performed a linear 
fit to this relation for each galaxy, estimating the blend 
center as the wavelength of minimum flux between 8180 
A and 8230 A. With respect to our best fit, we find a 
scatter in Nal of 0.08 A for NGC 1023 and 0.06 A for 
NGC 2974. These values are indeed comparable to our 
combined random and systematic errors in Nal, which 
range from 0.05-0.08 A for NGC 1023 and 0.04-0.07 A 
for NGC 2974. If anything, this test indicates that we 
have assessed our systematic errors conservatively. 

5. RESULTS: SPATIAL VARIATION IN LINE DEPTHS 

The radial variations in the Nal, Call, and FeH spec¬ 
tral features are illustrated in Figure |4] after convolv¬ 
ing binned spectra to the same rest frame and veloc¬ 
ity dispersion. In both galaxies, the Nal and Call fea¬ 
tures become visibly shallower toward larger radius. At 
a ~ 200kms _1 , the Nal doublet (8183 and 8195 A) is 
unresolved, and blended with a TiO band at 8205 A. 
In the left panels of Figure [4] it is apparent that as the 
Nal blend becomes shallower toward large radii, its cen¬ 


ter shifts toward redder wavelengths, an indication that 
the Nal feature is weakening more rapidly than TiO. 
This trend is qualitatively consistent with variations in 
sodium abundance or the stellar IMF; we discuss both 
possibilities in i |7.2l 

From Figure^] it is evident that the trend toward shal¬ 
lower FeH at large r is present in both the lin e region 
and the re d pseudo -continuum region defined bv lCvl) l~2l . 
While the[C vD12l definition highlights the deepest part 
of the FeH bandhead, their red pseudo-continuum still 
includes contributions from FeH and possibly TiO. To 
ensure that radial variations in the pseudo-continuum 
near FeH are not dominating our index measurements, 
we have tested three alternative variants of the FeH in¬ 
dex, including two that extend the line region to 9962 A. 
Details are provided in Appendix [Bj All variants yield 
similar trends in FeH with respect to r near the centers 
of NGC 1023 and NGC 2974. 

Radial trends in measured line indices are shown in 
Figure 0 for NGC 1023 and Figure l for NGC 2974. 
We group the H/3, bTiO, and TiOO.89 indices as indi¬ 
cators of age and temperature; the CNi, C24668, and 
(Fe) indices as indicators of C, N, or Fe abundance 
the Mg b, MgI0.88, and Ca4227 indices as indicators of 
a-process elements; and the Nal, Call, and FeH indices 
as IMF-sensitive indicators. Nonetheless, we stress that 
variations in the underlying stellar population have de¬ 
generate effects on multiple indices, and no single index 
or set of indices maps directly to a single stellar popula- 

6 We define (Fe) = 0.5(Fe5270 + Fe5335), followingI Trager et alj 
11200031) . 
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Figure 4. Comparison of spectra at multiple radial bins in each galaxy, featuring the gravity-sensitive absorption features Nal (left), Call 
(middle), and FeH (right). For visualization purposes only, we have divided each spectrum by a third-order polynomial. Each spectrum has 
been convolved to a velocity dispersion matching the central radial bin. Top row: NGC 1023, convolved to a — 230kms — 1 . Bottom row: 
NGC 2974, convolved to a = 245kms —1 . For FeH, light and dark rectangles indicate the line and pseudo-continuum regions, respectively. 
We also show a telluric emission spectrum (blue) and atmospheric transmission spectrum (grey), indicating the location of strong sky 
features. For each galaxy, the emission spectrum is scaled to the continuum level of the central bin, with additional scalings noted in panels 
c and f. The transmission spectrum is scaled by an additional factor of 0.3 in all panels. Both sky spectra match the rest frame of the 
central bin. 


tion property. In particular, we will inspect the meaning 
of Nal, Call, and FeH more carefully in 1 17.21 
All of the absorption features (except FeH in NGC 
2974) weaken toward large radii, as expected for galaxies 
harboring metallicity gradients. Yet there are notewor¬ 
thy differences between the rate of decline for different 
indices. The Nal and CNi features exhibit steep gradi¬ 
ents that appear nearly constant in log(r), while other 
features decline less steeply on average, and/or turn over 
toward a flat profile in the central kpc. Interestingly, 
the Wing-Ford FeH band does not mirror the steep ra¬ 
dial trend in Nal, although both are dwarf-sensitive fea¬ 
tures. In both galaxies, our measurements of FeH are 
consistent with uniform strength out to 0.2r e ff. Beyond 
this radius, NGC 1023 shows a gradual decline in FeH 
strength. NGC 2974 shows subtle evidence for increas¬ 
ing FeH strength toward large radii, but our outermost 
points for this galaxy are badly contaminated by telluric 
emission. The Ca4227, Mg b, and (Fe) indices all exhibit 
similar radial behavior, with a decline of 15%-25% per 
dex in r. Our trends in H/3, M g b, and (Fe) are broadly 
similar to those measured by iKuntschner et al.1 (|2006h 
with SAURON integral-field data out to r w 20". One 
exception is our measurement of increasing H/3 toward 
the center of NGC 1023. The absence of an H/3 gradient 


in the SAURON map may reflect differences in our re¬ 
spective methods for removing emission lines of ionized 
gas. Also, our measured Mg b values in both NGC 1023 
and NGC 297 4 vary more steepl y than the average trends 
displayed by IKuntschner et al.l (120061 1, though noise in 
their two-dimensional maps hinders a direct comparison. 

We display ratios of selected line indices in Figures 0 
and [HI For these figures, we have switched to a linear 
scale in radius so as to emphasize rapid changes within 
the central region of each galaxy. We have also removed 
the telluric absorption term from our error bars in Nal, 
as justified in 44.11 In Figure [7] we compare the Nal index 
to six other species. In every case except for CNi, the rel¬ 
ative strength of Nal increases toward the galaxy center, 
with a particularly steep rise in the innermost 0.1r e ff. 
This region turns ou t to be very si mil ar to the ry ,ff/8 
aperture size used by Ivan Dokkurn fc Conrovl (120121) . 

In Figure [5] we compare the radial variation of IMF- 
sensitive and a-element indices, relative to (Fe). Re¬ 
markably, the Nal and FeH indices show opposite trends 
with respect to (Fe) in the central 0.1 r e g. As in Fig- 
ureEJy NaI/(Fe) rises dramatically in the central 0.1r e ff 
of each galaxy. On the other hand, FeH/(Fe) decreases 
by ~ 10-20%. Even considering the large uncertainties in 
our measurements, the deviation between Nal strength 
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Figure 5. Selected line index strengths vs. radius in NGC 1023, corresponding to a = 230kms —1 . Vertical errors include systematic 
effects discussed in Appendix [A] and horizontal error bars indicate the radial bin sizes. Panels are ordered by wavelength and color-coded 
according to each index’s primary use as a diagnostic: grey (panels a, c, and g) for indicators of C, N, and Fe abundance; blue (panels 
b, f, and j) for a-element indicators; purple (panels d, e, and k) for stellar age and effective temperature indicators; and red (panels h, 
i, and 1) for IMF indicators. However, we stress that simultaneous modeling of multiple indices is necessary to quantitatively assess the 
contributions of star formation history, IMF, and abundances of individual elements. Most panels have units of A for equivalent width, 
except for CNi and bTiO (magnitudes), and TiOO.89 (ratio of blue to red pseudo-continuum). The y-axis in each panel scales from 0 to 1.3 
times the maximum line depth. Magenta diamonds in panel e represent H (3 measurements extracted without any correction for emission 
lines. The last data point in panel 1 (FeH, dashed error bar) is severely compromised by sky emission. Tabulated index values are available 
as online data. 
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Figure 6. Selected line index strengths vs. radius in NGC 2974, corresponding to a = 245kms —1 . Color codings and equivalent width 
units are the same as in Figure [5] The y-axis in each panel scales from 0 to 1.3 times the maximum line depth. The last two data points 
in panel 1 (FeH, dashed error bar) are severely compromised by sky emission. Tabulated index values are available as online data. 
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Figure 7. Ratio of Nal to other selected line indices as a function of radius in NGC 1023 (green) and NGC 2974 (blue). The spatial 
binning is the same as Figures 15151 but plotted on a linear scale. In panels a-e, the outermost spatial bins extend to 1.6r e ff for NGC 1023 
and 0.8r e ff for NGC 2974. For panel f (Nal/FeH), the outermost bins extend to 1.0r e ff for NGC 1023 and 0.26 r e ff for NGC 2974. Grey 
shaded areas indicate the scale of the single-aperture measurements bv lvan Dokkum fc Conrovl (|2012f) . corresponding to a physical radius 
of 320 pc in NGC 1023 and 480 pc in NGC 2974. Error bars include systematics from measuring kinematics and removing emission lines, 
but do not include a component for telluric correction of Nal (see Appendix I A. 3[) . For panel b, we measured the CNi equivalent width in 
A before comparing to Nal. 


and FeH strength appe ars t o be significant (Figure [7f). 
As discussed further in m this behavior can arise from 
a strong gradient in sodium abundance, whereas gradi¬ 
ents in the low-mass IMF slope would cause Nal and 
FeH to simultaneously increase or decrease. Although 
IMF variations still may be possible within these galax¬ 
ies, the opposing behavior of Nal and FeH constrains the 
magnitude and functional form of the IMF variations rel¬ 
ative to abundance variations. 

The near constancy of Mg b/(Fe) at all radii suggests 
that [a/Fe] is uniform or varies mildly with radius. The 
Mg b index is modestly sensitive to carbon, and our 
inferred gradients in [C/Fe] (' H7.3D ultimately allow for 
a shallow decrease in [a/Fe] toward large radii in each 
galaxy. In contrast to Mg b or Ca4227, Figure [S] shows 
an abrupt downturn in CaII/(Fe) interior to ~ 0.1r e ff. 
This steep central trend could arise from a lower fraction 


of giant stars in the very center of th e galaxy, or from an 
increase in sodium abundance (e.g., lCvD12f h Both ef¬ 
fects predict a simultaneous increase in Nal strength. In 
m we argue that sodium abundance drives the radial 
variations in Nal and Call. 

While the Nal gradients in NGC 1023 and NGC 2974 
are largely consistent, the CNi index shows mild discrep¬ 
ancies between the two galaxies. The overall gradient in 
CNi is approximately —0.06 mag per dex in r for NGC 
1023, versus -0.04 mag dex" 1 for NGC 2974. This dif¬ 
ference corresponds to a shallow decrease in Nal/CNi 
toward the center of NGC 1023, versus a shallow increase 
for NGC 2974 (Figure [7b) . 

6 . COMPARISON WITH SPS MODELS 

To complement our qualitative interpretation of the 
line index gradients in NGC 1023 and NGC 2974, we 
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Figure 8. Ratio of IMF-sensitive indices (Nal, Call, and FeH) to (Fe) as a function of radius in (a) NGC 1023 and (b) NGC 2974. We 
also illustrate the radial trends in Mg b/(Fe) and CNi/(Fe). The spatial binning is the same as Figures 15161 but plotted on a linear scale. 
T he outermost spatial bin i n pan el (a) extends to 1.6r e ff. The grey shaded area indicates the scale of the single-aperture measurements 
by Ivan Dokkum &; Conrovl (120120 . 


compare our measurements to model spectra bv lCvD12t 
Though rigorous fitting is reserved for future work, our 
comparison serves to highlight cases where strong radial 
variations in absorption line indices reveal underlying 
changes in abundance ratios or the stellar IMF. 

In Figure [9] we illustrate the radial variation of NGC 
1023 in index-index space, for H/3, Nal, Call, and FeH 
versus (Fe). For comparison, each panel includes vec¬ 
tors indicating the isolated effects of age, total metallicity 
[Z/H], [ct/Fe], and [C/Fe] on model spectra. Figure [TOl 
presents the same comparisons for NGC 2974. In each 
panel of Figures [9] and [TO] we display two “families” of 
model s extrap olated from publicly available model spec¬ 
tra by !CvD12l . In our construction, a family of models 
samples [Na/Fe] values of 0, +0.5, and +1.0, and four 
IMF variants: Chabrier, Salpeter (a = 2.35), and single 
(unimodal) power laws with a = 3.0 and a = 3.5. Each 
IMF in the lCvD12l models is integrated down to 0.1 M & . 
All parameters except for [Na/Fe] and the IMF are the 
same for a given family. For each galaxy, we construct 
one family to approximate the line indices in our central 
bin, and a second family at lower metallicity to approxi¬ 
mate the line indices at large radii. Instead of adjusting 
[Fe/H] in isolation, we find better overall agreement by 
varying [Z/H], a parameter available in more in recent 
models (Charlie Conroy; private communication). 

In order to test an appropriate range of ages, [Z/H], 
and other abundances, we construct a large number of 
model spectra ba sed on e xtrapolations from the baseline 
grid of models by [Cv ml Starting from a given age and 
IMF, the abundances are applied as multiplicative re¬ 
sponse functions, originally modeled for a Chabrier IMF 
at 13.5 CyiQ- We extrapolate response functions from 

' The response function for [Z/H] in the more recent models is 


the following baseline parameters supplied by [Cy DU 
[Z/H] G { -0.3, 0, +0.3 }; [a/Fe] G { 0, +0.2 }; [C/Fe] G 
{ -0.15, 0, +0.15 }; [N/Fe] G { -0.3, 0, +0.3 }; [Na/Fe G 
{ -0.3, 0, +0.3 }; and [Ca/Fe] G { -0.15, 0, +0.15 }. We 
then convolve each extrapolated spectrum from the na¬ 
tive resolution of the models to 230 kms -1 or 245 kms -1 , 
for comparison to data from NGC 1023 or NGC 2974, re¬ 
spectively. We compute line indices for each model spec¬ 
trum using the same procedure as for the galaxy spectra. 

To estimate the stellar population parameters corre¬ 
sponding to a particular galaxy spectrum, we compare 
the indices H/3, (Fe), Mg b, CNi, C24668, and Ca4227 to 
our set of extrapolated models. These indices are rela¬ 
tively insensitive to variations in [Na/Fe] or the IMF. We 
first find an age and total metallicity that approximately 
matches H/3 and (Fe) and then iteratively adjust the age 
and all abundance parameters above (except [Na/Fe]) 
until the resulting model spectra nearly match the values 
of all six indices. We emphasize that our method is not 
designed to produce a statistically rigorous fit, but rather 
to select families of models approximating the galaxy at 
small and large radii, so that we may visualize the effects 
of [Na/Fe] and the IMF on the near-infrared gravity- 
sensitive features. In the following section we examine 
the inferred abundances and possible IMF variations. 

7. DISCUSSION 

Variations in stellar masses, ages, and abundance ra¬ 
tios impose degenerate effects upon individual line in¬ 
dices. Inferring these physical properties demands an 
intricate comparison between observed data and stel¬ 
lar population and stellar evolution models. Herein we 
have attempted to present our measurements with suf- 

based on a Kroupa IMF. 
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Figure 9. Index-index trends in NGC 1023 and in SPS models, comparing different indices to (Fe). The thick colored line in each 
panel traces the radial variations in NGC 1023, with symbols marking our measurements for each radial bin. Overplotted are two model 
grids from ICvT*) I 'A - In each grid we vary the IMF (symbol size; solid line) and [Na/Fe] (symbol shape; dashed line). Close d grey s ymbols 
represent models at 10.5 Gyr, [Z/H] = +0.05, [a/Fe] = +0.3, [Ca/Fe] = +0.2, [C/Fe] = +0.4, and [N/Fe] = +0.15. The ICvD12i models 
with these parameters have similar Lick index values to our central bin for NGC 1023. Open red symbols represent models with [Z/H] 
= -0.5, corresponding to our outermost bins for NGC 1023. For simplicity, we have not varied age, [a/Fe], [Ca/Fe], [C/Fe], or [N/Fe] 
between the grey and red model grids. The labeled vectors represent the effects of varying age from 10.5 Gyr to 13.5 Gyr, [Z/H] by +0.55 
dex, [a/Fe] by +0.2 dex, and [C/Fe] by +0.3 dex. Our adjustments for total metallicity [Z/H] use recently updated SPS models (Charlie 
Conroy, private communication). 


ficient transparency to support future analyses employ¬ 
ing a wide range of modeling assumptions. With the 
caveat that rigorous interpretation requires careful mod¬ 
eling, we shall discuss some qualitative trends in the 
relative strengths of different line indices as a function 
of radius, in light of previously established connections 
to physical stellar prop ertie s. We discuss gradients in 
age and metallicity in Ef7.ll and then examine possible 
origins of the steep variations in Nal ( H7.2\) and CN/ 
(D, particularly in the context of star-forming pro¬ 
genitors of early-type galaxies like NGC 1023 and NGC 
2974 (1 17.41) . In 47.51 we examine each galaxy for photo¬ 
metric or kinematic signposts near 0.1 r e g, where the line 
indices and inferred stellar populations exhibit a sharp 
transition. Finally, in 47.61 we compare our findings to 
other recent investigations of radial IMF variations in 


early-type galaxies. 

7.1. Age and Metallicity Gradients 

Radial variations in H/3 suggest age gradients in both 
galaxies, although the H/3 index is sensitive to abundance 
variations as well as age (Figure [9^)- Based on the pro¬ 
cess outlined in ([SI we infer that NGC 1023 and NGC 
2974 both have old stellar populations at large radii, 
matching 13.5 Gyr models. Their centers are slightly 
younger: « 10.5 Gyr for NGC 1023 and ~ 11 Gyr for 
NGC 2974. As illustrated by the vectors in Figures [9[ 
age variations in this range have a weaker impact 
than [Z/H] or [Na/Fe] for all indices except H/3. We 
note that our ages are derived after removing emission 
lines from the galaxy spectra. The presence of spatially 
extended emission lines alongside an old stellar popu- 
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Figure 10. Index-index trends in NGC 2974 and in SPS models, comparing different indices to (Fe). The thick colored line in each 
panel trac es the radial variations in NGC 2974, with symbols marking our measurements for each radial bin. Overplotted are two model 
grids from [UvD121 . In each grid we vary the IMF (symbol size; solid line) and [Na/Fe] (symbol shape; dashed line). Closed g rey symbols 
represent models at 11.0 Gyr, [Z/H] = -0.1, [a/ Fe] = +0.3, [Ca/Fe] = +0.1, [C/Fe] = +0.3, and [N/Fe] = +0.2. The lCvD12l models with 
these parameters have similar Lick index values to our central bin for NGC 2974. Open red symbols represent models with [Z/H] = -0.6, 
corresponding to our outermost bins for NGC 2974. For simplicity, we have not varied age, [cn/Fe], [Ca/Fe], [C/Fe], or [N/Fe] between the 
grey and red model grids. The labeled vectors represent the effects of varying age from 11.0 Gyr to 13.5 Gyr, [Z/H] by +0.5 dex, [a/Fe] 
by +0.2 dex, and [C/Fe] by +0.2 dex. 


lation is consistent with photoionization from post-AGB 
stars, with smaller contributions from extre me horizontal 
branch stars and low-m ass X-ray binaries ( jBinette et al.1 
[1991 ISarzi eFalllmft . 

Both galaxies exhibit strong gradients in total metal- 
licity. In NGC 1023 [Z/H] declines from ss +0.05 in 
the central bin to ss -0.5 at r ~ r e ff. In NGC 2974 we 
measure [Z/H] ss -0.1 at the center and [Z/H] « -0.6 
at r ~ 0.6r e ff. The variations in [a/Fe] are shallower, 
from central values « +0.3 to outer values « +0.1 in 
both galaxies. We find that [Ca/Fe] is slightly less en¬ 
hanced than other a-el ements, similar to trends reported 
in other systems (e.g.. iThomas e t ah 20031: IGraves et all 

120071 : IWorthev et al.1 201 lib Examination of the CNi and 

C24668 index strengths suggests gradients in [C/Fe] from 
~ +0.1 at large radii to ~ +0.4 and ~ +0.3 at the centers 


of NGC 1023 and NGC 2974, respectively. In contrast, 
we infer [N/Fe] « +0.2 at all radii. Abundances of car¬ 
bon and nitrogen are discussed further in < 17.31 

We caution that comparing our inferred [Z/H] values 
for NGC 1023 and NGC 2974 to other measurements re¬ 
lies on a consistent definition of total metallicity. In the 
models we have employed, [Z/H] and other abundances 
are adjusted as independent variables, whereas in other 
cases [Z/H] may be a secondary quantity estimated from 
abundances such as [Fe/H] and [Mg/Fe]. In particular, 
secondary estimates of [Z/H] and [a/Fe] depend heav¬ 
ily upon oxygen abundance, which is difficul t to infer 
direc tly from stellar absorption features (e.g., iSchiavonl 
[ 2001 / 

Having approximated the ages and abundances of our 
inner and outer bins for NGC 1023 and NGC 2974, we 
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Figure 11. Index-index trends in NGC 1023 and in SPS models, comparing IMF-sensitive indices. Data, model grids, and vectors for 
individual parameter variations are defined as in Figure [9] 
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Figure 12. Index-index trends in NGC 2974 and in SPS mode ls, comparing different IMF-sensitive indices. Data, model grids, and 
vectors for individual parameter variations are defined as in Figure [Tril The age vector is not displayed, as the difference between 11.0 Gyr 
and 13.5 Gyr is even smaller than the [C/Fe] vector for these indices. 


can select appropriate families of models to isolate the ef¬ 
fects of [Na/Fe] and IMF variation on various absorption 
indices. In spite of the weak trends in age, [ct/Fe], [C/Fe], 
and [Ca/Fe] reported above, we choose in Figures IfllfTill 
to illustrate grids that differ only in [Z/H]. The impacts 
of other abundance patterns are smaller in magnitude 
and are instead represented by the vectors in each figure 
panel. For both grids the adopted age, [ct/Fe], [C/Fe], 
[N/Fe], and [Ca/Fe] match the approximate values for 
the galaxy center. With these abundances supplying a 
baseline, each grid then represents a family of models 
with varying [Na/Fe] and IMF. 


We now consider two absorption indices that depart 
from the trend of gradual radial decline: Nal and CNi. 
As illustrated in Figure [3 the Nal index exhibits a much 
steeper gradient than other features in the central few 
hundred pc of each galaxy. Figures fill and fl2l also indi¬ 
cate a clear break in the behavior of Nal for the central 
3-4 bins in each galaxy, relative to the trends in Call 
and FeH. As discussed further in m this break appears 
consistent with an abrupt rise in sodium abundance in¬ 
terior to ~ 0.1r e ff, and limits radial IMF variations to 
stellar masses above ~ 0.4 M 0 . Surprisingly, CNj is the 
only feature showing similar radial variation to Nal (Fig- 
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ure Eb). We have considered both carbon and nitrogen 
abundance variations as a possible explanation for the 
steep gradient in CNi and discuss these possibilities in 

7 . 2 . IMF vs. Sodium Abundance Variations 

The Nal feature is primarily sensitive to surface tem¬ 
perature, surface gravity, and sodium abundance. The 
second effect makes it a strongly dwarf-sensitive feature 
in uniformly old stellar populations. Thus it is tempt¬ 
ing to interpret our observations as a steep gradient from 
a bottom-heavy IMF at the galaxy center to a shallower 
IMF slope at r A, 0.1 r e g (~ 300 pc). However, this inter¬ 
pretation must be reconciled with the relatively mild de¬ 
cline in FeH and the opposing behavior of Nal/(Fe) and 
FeH/(Fe). While the Nal and FeH features both peak in 
sensitivity for stellar masses m < 0.2 Mq, the sensitivity 
of FeH declines much faster toward higher masses, such 
that Nal is more than twice as sensitive to the num - 
ber of stars with m > 0.4 Mq (see Figure 17 of lCvD12l) . 
Thus a strong radial trend in Nal but not FeH could 
reflect an IMF gradient whose “bottom-heavy” nature is 
only expressed above ~ 0.4 Mq. A bimodal IMF whose 
slope only varies above ~ 0.5 M 0 has been shown to 
agree with the dynamical masses of early-type galaxie s 
(e.g.. lLa Barbera et al.ll20l3ll2015HSpiniello et al.ll20T4l) . 
Yet the magnitude of the variation in Nal would de¬ 
mand drastic IMF variation on sub-kpc scales if the IMF 
slope above ~ 0.4 M 0 were the sole driver. The bTiO in¬ 
dex also varies less steeply than Nal, although its strong 
temperature-sensitivity and the overlapping Mg4780 line 
leave some doubt over whethe r it is a valid IMF indi- 
cator (e.g.. iServen et al.l 120051 lLa Barbera et al.l 120131 
I SomielloetllJ l 201 4lb - 

Varying sodium abundance offers an alternative expla¬ 
nation for the steep variation of the Nal feature inside 
~ 0.1 r e ff. Indeed, the outer regions of both galaxies align 
with nearly solar [Na/Fe] in Figures [Til and fl2l while the 
central index strengths indicate [Na/Fe] between +0.5 
and +1.0. Due to its prominent role as an electron donor, 
sodium impacts other features via the atmospheric elec¬ 
tron pressure. In particular, an increase in [Na/F e] will 
drive a mild decrease in ionized calcium (tCv Dial . This 
is qualitatively consistent with the flattening we observe 
in the Ca II index near the center of NGC 1023 and NGC 
2974, and the corresponding decline in CaII/(Fe). 

Comparing Nal and FeH in Figures UTb and fl2b. we 
find that both galaxy centers occupy the same locus on 
the central model grid: [Na/Fe] is just below +1.0 and 
the unimodal IMF slope lies between 2.35 and 3.0. The 
trends in Nal versus Call are more difficult to interpret, 
as IMF and [Na/Fe] variations are more degenerate for 
this pair of indices. Furthermore, the Call predictions 
from the SPS models are sensitive to our estimate of 
[Ca/Fe]. Allowing for an error ~ 0.5 A in the horizontal 
placement of the model grids in Figures [TTh and I12h , 
we again find that both galaxy centers are consistent 
with a sodium enhancement of approximately 1.0 dex 
and an IMF slope between 2.35 and 3.0. At large radii, 
the Nal and Call strengths of both galaxies (as well as 
FeH in NGC 1023) are roughly consistent with a Chabrier 
IMF. We note that while t he Cha brier form is directly 
employed in the models by !CvDl2 . it is not a unimodal 
power law, and the apparent transition from a > 2.3 to 


Chabrier may not fully reflect the unimodal or bimodal 
form of the IMF at different radii. 

Spatially resolved measurements of the Na D feature 
near 5890 A would provide additional evidence for or 
against gradients in [Na/Fe]. Unfortunately, our settings 
for LRIS did not c over Na D in NGC 1023 or NGC 2974. 
iJeong et al.l (12013H constructed a stacked SDSS spectrum 
of early-type galaxies with strong Na D absorption and 
reported that the Na D line strength requires super¬ 
solar [Na/Fe] as well as [Z/H], with IMF slope having 
a relatively minor effect on the Na D feature. Although 
their analysis provides circumstantial evidence in favor of 
[Na/Fe] driving trends in Nal, their measurements leave 
room for IMF variations in a subset of their galaxies, 
which are stacked over cr ~ 150-300 km s -1 . 

Sodium is produced primarily through the Ne-Na 
chain, which can be acti vated during core burn¬ 
ing in massive star s (e.g.. IWooslev fe Weaved 119951 : 
iDecressin et al.l l2007t iKobavashi et al.l I2006L I2011D or 
at the base of the convective envelope in intermediate- 
mass stars on the asymptotic giant branch (AGB), 
a process known as hot bottom burning (e.g, 
Cottrell^^a^^CostaLJl^Sl : iDenise nkov &; Dcniscnkova 

1990HVentura fe P’Antona l2008al lbl: lKarakadl2010l) . Al¬ 

though the dredge-up of heavy elements to the outer en¬ 
velope of AGB stars is most efficient at low metallicities, 
solar-metallicity AGB stars may still pro duce and eject 
non-negligible quan t ities of sodiu m (e.g.. lMowlavilll999l : 
iKarakas et al.l[20021 : iHerwigl1200511 . 

In addition to IMF and [Na/Fe] variations, [CyDT2| de¬ 
termined that Nal could be strengthened by decreasing 
the number of horizontal branch and AGB stars, or the 
number of extremely cool giants (M7III). Yet both of 
these effects substantially weaken the TiOO.89 feature, 
which we observe to be nearly constant over the radii 
where Nal varies most steeply. While our interpreta¬ 
tion of Nal and other features is self-consistent within 
the framework of our adopted SPS models, some recent 
studies have cautio ned against usi n g sodiu m indices to 
investigate the IMF. lSpiniello et al . ( t2015a l noted incon¬ 
sistencies between pred icted sodium depths from differ¬ 
ent SPS models, while iSmith et al.1 (l2015blf found dis¬ 
agreement between the total mass-to-light ratios in two 
lensing galaxies and the IMF and [Na/Fe] values inferred 
from the Na D and 1.14 ^m Na features. 

7 . 3 . Carbon vs. Nitrogen Abundance Variations 

The CNi index is sensitive to carbon and nitrogen 
abundance and exhibits a strong radial trend in both 
galaxies. Our other carbon-sensitive feature, C24668, 
does not decline as steeply as CNi from r = 0 to 
r ~ 0.1r e ff, and thus it is tempting to interpret the 
trend in CNi as a decrease in nitrogen abunda nce. Yet 
our comparison with the SPS models bv lCvD12l suggests 
that a gradient in [C/Fe] is better able to reproduce the 
observed trends in CNi and C24668. This counterintu¬ 
itive result occurs in part because CNi is more sensi¬ 
tive to [Z/H], such that the metallicity gradient in each 
galaxy drives a seemingly exaggerated gradient in CNi 
strength. Nonetheless, the absolute values of CNi rela¬ 
tive to (Fe) indicate that nitrogen is mildly enhanced at 
all radii, with [N/Fe] ~ +0.2. Neither CNi nor C24668 
are sensitive to sodium abundance or IMF variations. 

Carbon is produced by helium fusion in the cores of 
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stars and released into the interstellar medium (ISM) 
via stellar winds or Type II supernovae. Models of 
massive stars (m > 10 Mg) produce sizable carbon 
yields at solar and slightly lower metallicities (e.g., 


yields at solar and slightly lower metallic ities (e.g., 
Mevnct fc Maederjl2002bldDrav ct ahll2003l : [brav fc Toud 


20031 ). In intermediate-mass stars the final carbon yield 
is sensitive to initial mass and metallicity, as hot bot¬ 
tom burning efficiently converts carbon to nitrogen. 
At solar metallicity, stars with m ps 2.5-4 Mq trans¬ 
port carbon-rich material to the surface during third 
dredge-up events, without reachin g sufficient tempera¬ 
tures for hot bottom burning (e.g. , IRenzini fc Volilll981fc 
Gavilan et al.1 120051: Herwia 120051 : iKarakas fc Lattanziol 
2014 IVentura et al l 20151) . Abundance patterns in the 


Milky Way disk and nearby star-forming galaxies ■ par¬ 
ticularly an upturn in [C/O] toward solar metallicity - 
have prompted differing interpretations of the dominant 
source of carbon enrichment in these environments. On 
one hand, the absolute quantities of carbon and oxygen 
are consistent with open-box models of chemical evo¬ 
lution in st ar-forming galaxies, em p loying massive sta r 
yields (e.g., iGustafsson et all 119991 IHenrv et all l2000h . 
Conversely, observations of nearly constant [C/Fe] over 
~ 1 dex in [Fe/H] suggest that the ISM becomes carbon- 
enriched on timescales similar to Type la supernovae, 
point ing to intermediate-mass stars as the main supp lier 
(e.g., iChiappini et al. I [20031 : iBensbv fc Feltzin3l2006l) . 

Nitrogen is produced via the CN cycle. Whereas pri¬ 
mary nitrogen production is driven by the mixture of 
hydrogen- and helium-burning regions in the first gen¬ 
eration of stars, the nitrogen abundance of moderate- 
and high-metallicity stellar populations is dominated by 
secondary production from later generations whose inte¬ 
riors are already seeded with carbon. Similar to sodium, 
nitrogen production may occur in the cores of massive 
stars or in th e stellar envelopes of AGB stars above 
A Mr* (e.g.. IRenzini fc V'olil 1 1 98Tt iWooslev fc Weaved 


1995HGavilan et ah 20051 Ventura fc D’Antonall2008al lbf 


Karakas fc Lattanziol 12014b Observations of [N/Fe] in 
the Galactic disk and halo - including globular clus¬ 
ters ranging from [Fe/H] ~ —2.5 to [Fe/H] ~ —0.5 - 
are consistent with contributions from both massive and 
intermediate-mass stars (e.g., Chiappini et al. 2005; 
Hirschi 2007; cf. Cohen et al. 2005). 

Models of fast-rotating massive stars can boost CNO 
yields by mixing the hydrogen shell and helium-burning 
core (e.g., iMevnet fc MaedeHl2002al Tbl). In contrast, hot 
bottom burning in AGB stars will drive a strong anti¬ 
correlation between [C/Fe] and [N/Fe] unless a second 
mechanism is responsible for producing one of these el¬ 
ements. Although this anti-correla tion has been ob¬ 
served in some globular clus ters le.g.. lCohen et al.ll2005l : 
IVentura fc D’Anton a 2008b), the radial variation of CNi 
and C24668 in NGC 1023 and NGC 2974 is incon¬ 
sistent with opposing gradients in [C/Fe] and [N/Fe]. 
Rather, the trend of rising [C/Fe] at constant [N/Fe] 
suggests separate sources o f carbon and nitrogen pro¬ 
duction. IHenrv et~aT. (200 0) reached a similar conclusion 
after compiling data from individual stars in the Galac¬ 
tic disk and halo, plus Galactic and extragalactic HII 
regions. 


7.4. Abundance Ratios in Composite Stellar Populations 


The stellar populations of early-type galaxies are typ¬ 
ically old, with high metallicities resulting from multiple 
generations of star formation and ISM enrichment in a 
deep potential well. Their high [a/Fe] ratios could arise 
from short star formation timescales or a top-heavy IMF, 
both of which rapidly seed the ISM with q-process ele¬ 
ments via Type II superno vae ('e.g.. iWorthev et ahlll992l : 
iThomas et al.l[l999l I2005D . Here we briefly explore sce¬ 
narios that could yield an excess of sodium and carbon in 
the centers of early type galaxies, with steep abundance 
gradients to larger radii. 

The most straightforward hypothesis for excess sodium 
and carbon is a larger fraction of the stars responsible for 
producing these elements. This could be the direct re¬ 
sult of an IMF with a shallower slope above ~ 10 Mq 
if sodium and carbon are produced in massive stellar 
cores. The form of the IMF at high densities is an on¬ 
going challenge for models of star formation, and may 
be especially sensitive to the initial gas density struc¬ 
ture and the role of turbul e nce in driving fr agmentation 
('e.g.. iChabrier et all 12014 iKrumhoIzlI2014D . However, 
IMF variations cannot be invoked to explain correlated 
gradients in [C/Fe] and [Na/Fe] if AGB stars are the 
dominant source of both sodium and carbon. This is be¬ 
cause the IMF above ~ 2 Mq sets the number of stars 
that experience hot bottom burning, which exerts oppo¬ 
site influences on carbon and sodium yields. 

Another factor in the observed excesses is the amount 
of enriched stellar ejecta accreted onto existing and newly 
forming stars, during the peak of star formation. For 
instance, the high stellar densities in globular clusters 
allow for prolific accretion onto existing stars, as long 
as the cluster is massive enough to retain low-velocity 
ejecta (e.g., Renzini 2008; Conroy 2012; cf. Fenner 
et al. 2004) Usin g the deprojection procedure of 
iCebhardt et all (11996D . we find that the luminosity densi¬ 
ties of NGC 1023 and NGC 2974 rise from 2-6 Lq pc~ 3 at 
0.1 r e ff to 200-400 Lq pc -3 at 0.01 r e g, comparable to the 
average densities of some globular clusters. An extremely 
high value of [Na/Fe] ps 1 dex has been measu red in 
the central 15 pc o f M31 dConrov fc van Dokkuiull2012bl : 
IZieleniewski et ahll2015D . while red giants in the Galactic 
bulge and Galactic globular cluster s both exhibit [Na/Fe] 
between 0 dex and +0-5 dex (e.g., iLecureur et al.1120071 : 
iRoediger et al.112014 I.Tohnson et all 12014 120151b Al¬ 
though support for this mechanism in globular clusters 
is still contentious, the rapid increase in sodium abun¬ 
dance toward these galaxies’ densest regions is at least a 
plausible consequence of pollution by stellar ejecta. 

Compiled measurements of sodium, carbon, and ni¬ 
trogen in other early-type galaxies do little to clarify 
the observations reported herein. Spectra of the cen¬ 
tral rs 0.6-1 kpc in early-type galaxies exhibit stronger 
CNi features and steeper variation of Na D versus (Fe) 
than composite stellar population models for e llipti¬ 
cal galaxies or the Galactic bulge (iTrager et all [2004 
iServcn fc Worthevll2010HTang et al.ll2014l) . iGraves et al.1 

l2007ir and lJohansson et al.l (120121) report slightly higher 
carbon enhancement than nitrogen en hancement in 
stacked spe ctra of SDSS galaxies, while iGreene et al.l 
(120131 [20151 ) report higher [N/Fe] than [C/Fe] for a com¬ 
bined sample of 82 galaxies. Although stacked early-type 
galaxy spectra allegedly range from 0 dex to +1 dex in 
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Figure 13. Major-axis photometry of NGC 1023. (a) Surface 
brightness pr ofile (dashed grey line), combining HST-based mea¬ 
sur ements bvJLauer et aLl (120051) and ground-based measurements 
by IScott et al.l (120091) . The solid black line illustrates the corre¬ 
sponding de projected luminosit y density, (b) V-I color and ellip- 
ticity e from fLauer et ahl <120051 ), overlaid with our measurements 
of NaI/(Fe). 


Figure 14. Major-axis photometry of NGC 2974. (a) Surface 
brightness pr ofile (dashed grey line), combining HST-based mea- 
su rements b y Laucr et al. (2005) and ground-based measurements 
bv IKrainovic et al.l (120051) . The solid black line illustrates the cor¬ 
responding d eprojected luminos ity density, (b) V-I color and ellip- 
ticity e from fLauer et ahl (120051 ), overlaid with our measurements 
of NaI/(Fe). 


[N/Fe], the variations between different studies are larger 
than the inferred trends within any individual sample, 
suggesting systematic errors in at least some measure¬ 
ments of [N/Fe]. 

7.5. Is 0.1r e ff a Special Radius? 

Even while the causes of excess sodium, carbon, and 
nitrogen are murky, we have shown that the most ex¬ 
treme stellar populations in NGC f023 and NGC 2974 
reside in their innermost regions. In both objects, the 
strengths of Nal, Call, and CNi relative to other in¬ 
dices change abruptly near 0.1 r e g, leading us to question 
whether this scale marks a transition in the galaxy’s pre¬ 
vious star formation environment, or a structural land¬ 
mark from the assembly of distinct progenitors. Seeking 
supporting evidence for unique behavior at 0.1r e ff, we 
have exa mined Hubble Space Telescope (HST) photom¬ 
etry from lLauer et al.l (l2005f) to examine each galaxy’s 
surface brightness profile, V — I color, and ellipticity out 
to 0.2r c ff. We also exami ned the more ex t ended sur- 
face brightness pr ofiles from IKrainovic et al.l (|2005fl and 
IScott et ahl (I2009H . as well as our measurements of v and 
<t for each radial bin. Wc illustrate the photometric data 
in Figures [Hi] and 1141 and the kinematic data in Ap¬ 
pendix IA.1I (Figure [15]) • 

Generally speaking, neither galaxy exhibits an abrupt 


feature near 0.1 r e fi in its stellar kinematics or broadband 
light. In NGC 2974, a is the dominant kinematic compo¬ 
nent inside 0.1r e fi, whereas rotational v dominates be¬ 
yond 0.2r e ff. NGC 2974 also has a central dust feature, 
which intersects the major axis near 0.05 r e ff. However, 
the absence of similar features in NGC 1023 suggests that 
they are not closely connected to the sharp abundance 
variations. 

7.6. Other Reports of IMF-Sensitive Index Gradients 

To date, few studies have sought to measure IMF gradi- 
ents within individual ea rly-type galaxies. Most recently, 
lLa Barbera et all (120151 ) fit radial trends in TiO0.89 and 
four temperature- and gravity-sensitive TiO features in 
a massive early-type galaxy (a ~ 300kms _1 ) at redshift 
z ss 0.057. Their measurements are consistent with a 
bimodal IMF whose slope above ~ 0.5 Mq declines from 
a centr al value a ~ 3 to a Milk y Way-like IMF beyond 
0.5 r e ff. El Barbera et al.l ( 2015 ) have measured constant 
FeH strength out to « 0.2 r e g (~ 1.6 kpc), which dis¬ 
favors a varying unimodal IMF. Our measurements of 
FeH in NGC 1023 and NGC 2974 exhibit a similar trend 
at r < 0.2r e ff, and probe physical scales several times 
smaller than the centra l bin of ~ 700 pc adopted by 
lLa Barbera et al.l (|2015l) . 
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IZiclcniewski et all (j2015l) have measured a rapid in¬ 
crease in Nal over the central 20 pc of M31, with no 
corresponding increase in FeH. Similar to our interpre¬ 
tation above, they infer a steep gradient in [Na/Fe]. The 
central 40 pc of M32 do not exhibit strong radial trends 
in Nal or FeH. 

IMartin-Navarro et al.l (l2015aT) have measured gradi¬ 
ents in IMF-sensitive features in three galaxies, using 
long-slit data covering 4500-10,000 A. Near the centers of 
their two high-mass galaxies, NGC 4552 and NGC 5557, 
they find that NalsDSS varies steeply relative to the to¬ 
tal metallicity indicator [MgFe]', similar to the trends 
we find above. Their lower-mass object, NGC 4387, 
shows very little radial variation in Nalgoss or i n inferred 
stellar population properties, aside from total metallic¬ 
ity. Whereas we have carefully accounted for telluric 
absorption near 819 0 A (see H3.21 and Appendix IA.3I) , 
IMartin-Navarro et all ( 2015a ) place less confidence in 
their telluric correction and exclude the NalsDSS index 
from their SPS models. They do not measure the FeH 
feature, and telluric emission restricts their assessment of 
the Call feature to the second line in the triplet. How¬ 
ever, they observe and model other features between 5800 
A and 6400 A that we do not access with LRIS: most 
notabl y the Na D index and tw o TiO indices. 

IMartin-Navarro et al.l (I2015af) have reported a strong 
IMF gradient in NGC 4552, resolved to ~ 0.1 r e s or 
~ 300 pc. They a dopt the same bimodal IMF form as 
lLa Barbera et all (120151) and likewise measure a Milky 
Way-like IMF at 0.7r e ff, while at the center of NGC 
4552 they find an extreme slope a « 4 above 0.5 Mq. 
They illustrate that NalsDSS and Na D exhibit similar 
gradients in NGC 4552 despite different sensitivities to 
[Na/Fe], and assert that IMF gradients mainly drive the 
variations in Nal, while Na D permits [Na/Fe] variations 
up to 0.25 dex. While the combined coverage of Na D 
a nd Nal provides use ful l everage , the modeling approach 
of IMartin-Navarro et all (I2015af) warrant s caution: they 
fit sele cted line indices to SPS models bv IVazdekis et al.l 
• 2012 :. after applying correction factors to adjust the 
measured indices from inferred [a/Fe] to solar abun¬ 
dances. The overt dependence on [a/Fe] is particularly 
troubling for NGC 4552, where the [a/ Fe] gradient re¬ 
ported bv IMartin-Navarro et al.l (f2015al ) is much steeper 
than the trends typically observed for early-type galaxie s 
(e.g.. iSanchez-Blazquez et ahlf2007l : iSnolaor et al.l ; 2010 ). 

Using the the same modeling approach and bimodal 
IMF form, IMartin-Navarro et al.l (l2015bl) reported a 
high-mass slope a « 4.0 out to > 1.5r e ff in the com¬ 
pact, high-er galaxy NGC 1277, with a mild trend toward 
a ss 3.5 in the central 0 .5 r B ff (600 pc). For this object, 
IMartin-Navarro et all ( 2015a ) found strong gradients in 
NalsDSS) Na D, and metallicity, compared to relatively 
weak trends in TiOi, TiC> 2 , and [Mg/Fe]. The results 
for NGC 1277 agree qualitatively with our observations 
of NGC 1023 and NGC 2974 on similar radial scales, 
and with widespread trends in [Z/H] and [a/Fe] (e.g., 
Tamura et al.l 1200(1 : iWeiimans et all 120091: iGreene et all 
20131I2015D . Finallv. lMartm-Navarro et al.1 (l2015d) have 
presented data, from 24 galaxies in the CALIFA survey 
(|Sanchez et al.|[2012D . binning each galaxy into several 
elliptical annuli. Although they find a significant cor¬ 
relation between the inferred IMF slope and metallic¬ 


ity of each spectrum, they do not compare any mea¬ 
surements directly with r. We therefore cannot assess 
whether galaxies in the CALIFA sample contain IMF or 
abundance gradients at the spatial scales we have probed 
f or NGC 1023 and NGC 2 974. 

lLa Barbera et al .i (l2013i ) assessed stacked spectra from 
SDSS and noted that the individual galaxy spectra en¬ 
closed varying fractions of r e g. They subdivided two 
of their stacked spectra (a = 100 kms -1 and a = 
200 km s -1 ) into narrow bins of r/ r e s and found minimal 
variations in the strengths of TiO indices and NalsDSS 
from aperture sizes of 0.3r e ff to 1.4r e ff. In other words, 
their a = 200 kms -1 stack exhibited stronger IMF- 
sensitive features at all enclosed radii. Still, we note that 
any radial gradients in abundances or the IMF would 
be diluted by the luminosity-weighted SDSS apertures. 
Our finding that the steepest gradients in NGC 1023 and 
NGC 2974 occur well inside 0.3 r e g m ay also bear upon 
the ab sence of gradients reported by ILa Barbera et al.l 

mm- 

Further advances in stellar template libraries and at¬ 
mospheric models will yield improved predictions for 
the integrated-light signatures of simultaneous variations 
in abundances and the IMF. In the meantime, models 
whereby abundance variations are restricted to [Z/H] or 
[a/Fe] must be employed with caution, in light of abun¬ 
dance gradients for indivi dual elements. An inter esting 
case study is provided bv lLa Barbera et all (120151 ). who 
attempt to rescale their observed index strengths to so¬ 
lar abundances, based on relations between a given in¬ 
dex strength, [a/Fe], and [Z/H]. Their approach is in¬ 
formed in part by their assessment of radial trends in 
[q/Fe] and [C/Fe ] . We speculate that the approach of 
lLa Barbera et al.l (j2015l) is better suited to the gravity- 
sensitive TiO features they measure than to the Nal or 
Call features discussed herein, since [Ca/Fe] and [Na/Fe] 
are known to deviate from trends in a/Fe]. We look for¬ 
ward to further work that may support this example and 
identify other regimes where IMF gradients can be ro¬ 
bustly disentangled from nuanced abundance variations. 

8. CONCLUSION 

We have used Keck/LRIS to analyze optical and near- 
infrared stellar absorption features along the major axis 
of two early-type galaxies with central a of 217 kms -1 
and 247 kms -1 . We have measured 13 line indices for 
species of H, C 2 , CN, Na, Mg, Ca, TiO, Fe, and FeH, 
from spatially resolved spectra covering 3100-5560 A and 
7500-10,800 A. We have examined each index on scales 
~ 100 pc (1") near the center of each galaxy, and in 
larger bins extending to 4.0 kpc (75" or 1.6r e ff) for NGC 
1023 and 3.0 kpc (30" or 0.8r eff ) for NGC 2974. 

In both galaxies, radial declines in the (Fe) index and 
multiple indices of Mg, Ca, and C suggest an overall 
metallicity gradient and shallower gradients in [a/Fe], 
[Ca/Fe], and [C/Fe]. However, the Nal index at 8190 
A exhibits significantly steeper gradients, particularly 
at r < 0.1 r e ff, or r < 300 pc. The FeH index at 9915 A 
mirrors the gradual decline of (Fe) rather than the steep 
decline in Nal. The data presented herein are among 
the first to track FeH as a function of radius, and to 
demonstrate different radial trends in Nal and FeH, even 
while both indices are sensitive to cool dwarf stars. 
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We interpret the steep gradient in the Nal index as re¬ 
flecting a rapid decline in [Na/Fe] over the central ~ 300 
pc of each galaxy. On similar scales, the Call index 
declines relative to (Fe) and other Mg and Ca indices, 
to a degree that would require a very strong gradient 
in [Na/Fe] if our interpretation is to match models by 
lCvD12l . IMF gradients may contribute to the respective 
trends of Nal, Call, and FeH in NGC 1023 and NGC 
2974, but only if the IMF is bimodal, with a similar slope 
for m < 0.4 Mg at all radii. The CNi index also varies 
rapidly near the center of each galaxy, and a qualitative 
comparison to stellar population synthesis models indi¬ 
cates that this is due to a large-scale gradient in [C/Fe], 
with uniform enhancement in [N/Fe]. 

Our study poses a number of outstanding issues to be 
pursued as future work. As we have emphasized above, 
the physical properties of stellar populations are highly 
degenerate with individual line indices, and some of the 
trends we have presented - most notably the relative 
variations in Nal, FeH, and (Fe) - defy simple qualita¬ 
tive arguments. Rigorous stellar population synthesis 
modeling will expand upon the qualitative comparisons 
we have performed herein and allow for more robust 
interpretations of the physical trends in NGC 1023 
and NGC 2974. Moreover, our search for radial IMF 
trends within two galaxies cannot fully inform claims of 
IMF variation over a large range of integrated galaxy 
properties. To this end we have observed five additional 
galaxies spanning a ~ 140-400 km s -1 , approaching 
the range explored by survey-driven investigations of 
IMF variation. These objects will also strengthen our 
understanding of internal trends near 0.1r e ff and clarify 
whether this scale truly marks a transition point for star 
formation physics in galaxies. 
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APPENDIX 

A: SYSTEMATICS IN LINE INDEX MEASUREMENTS 

Herein we derive our total systematic error e sys , which 
we have included above in H4 11 Figures [3] and [S] above 
include all systematic terms, while Figures [30 include all 
systematic terms except for the telluric absorption error 
in Nal. 


A.l: Kinematic Corrections 

In order to make an unbiased comparison between 
absorption line strengths at different radii, each spec¬ 
trum within a single galaxy must be shifted to the same 
rest frame and convolved to the same velocity disper¬ 
sion. For each spectrum, we measure v and cr us¬ 
ing the temp l ate reconstruct ion procedure pPXF from 
iCappellari fe Emsellcml 1 2004 ]. We use a subset of tem¬ 
plate sp ectra from t he MILES library of empirical stellar 
spectra dSanchez-Blazauez et al.f2006T) , and perform a fit 
over the wavelength range 3650-5400 A. We simultane¬ 
ously fit for the strengths of emission lines, assuming that 
gas and stars exhibit the same kinematics. Spectra on 
the LRIS red arm are assumed to match the kinemat¬ 
ics derived from blue-arm data for the same spatial bins. 
We display the radial profile of v and a for each galaxy 
in Figure [Tsl 

At high S/N, uncertainties in v and cr are domi¬ 
nated by systematic errors, particularly from template 
libraries that do not fully reproduce a galaxy’s under¬ 
lying stellar population. For this work, we do not at¬ 
tempt to directly assess our systematic errors in kine¬ 
matic moments, but rather test the variation in line index 
strengths over a conservative range of assumed v and cr. 
Previous investigations using similar wavelength cover- 



r/r eff 

Figure 15. Radial velocity and velocity dispersion profiles for 
NGC 1023 and NGC 2974, as measured with pPXF. At each radius 
we have averaged |u| and cr from opposite sides of the galaxy. We 
define v = 0 at our central bin for each galaxy. Error bars represent 
a constant systematic error of lOkms -1 in v and 20kms —1 in cr. 
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Table 3 

Systematic Errors in Line Indices 


Index 

Error from t;o 

Error from ctq 

Telluric absorption error 
NGC 1023 / NGC 2974 

Error from emission lines 
NGC 1023 / NGC 2974 

Adopted systematic error 
NGC 1023 / NGC 2974 

CNi 

0.4% 

0.9% 


0.9 % / 1.4 % 

1.3 % / 1.7 % 

Ca4227 

1.3% 

1.9% 

- 

- 

2.3 % 

C 2 4668 

0.4% 

0.4% 

- 

- / 0.4 % 

0.6 % / 0.7 % 

bTiO 

4.1% 

1.4% 

- 

- / 5.2 % 

4.3 % / 6.8 % 

H/3 

1.4% 

0 .2% 

- 

6.0 % / 6.0 % 

6.2 % / 6.2 % 

Mg b 

0.4% 

0.9% 

- 

0.7 % / 1.0 % 

1.2 % / 1.4 % 

Fe5270 

0.3% 

1.4% 

- 

- 

1.4 % 

Fe5335 

0.7% 

2.7% 

- 

- 

2.8 % 

Nal 

0.025 A 

5.5% 

8.8-13.4% / 1.2-3.8% 

[0.025 2 

+ e 2 el +(0.055 EW N al) 2 ] 1/2 A 
= 11.9-23.1% / 7.6-11.8% 

Call 

1 .1% 

0.7% 

_ 

_ 

1.3 % 

MgI0.88 

2 .1% 

5.6% 

- 

- 

6.0 % 

TiOO.89 

0.03 % 

0.06 % 

- 

- 

0.07 % 

FeH 

2 .6% 

2.3% 

- 

- 

3.5 % 


Notes: In most cases the adopted systematic error is e sy s = (e?o + £ cto) 1,/2 ! where e„o and e CT o are the error terms from vq and oo. 
For indices impacted by telluric absorption or galaxy emission lines, the adopted systematic error is e sys = (eJo+Eo-o+Etei+^gas) 1 ' 2 , 
where e te i is the error resulting from uncertainties in the atmospheric transmission spectrum, and e gas is the error derived from 
different settings for emission line removal. For Nal, the adopted systematic error ranges from 11% of the equivalent width near 
the galaxy center to 16% of the equivalent width at large radii. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 0.6 

Nal EW (A) Nal EW (A) 


Figure 16. Top: Error in Nal feature strength for test spectra in NGC 1023 and NGC 2974, where each spectrum was shifted to rest 
wavelengths from a range of velocities vq. Uncertainty in the systemic velocity produces a nearly constant error over a large range of line 
depths. The dashed line corresponds to an average error of 0.025 A, from trials sampling a conservative distribution of Do- Bottom: error 
in Nal strength for test spectra convolved to 230kms^ 1 (NGC 1023) or 245 krn s 1 (NGC 2974) from a range of initial velocity dispersions 
rjQ. Uncertainty in the measured velocity dispersion produces an error that increases with line depth. The dashed line corresponds to an 
error of 5.5%. 


age and a variety of stellar templates have fo und total er¬ 
rors 10 kms -1 i n v and ($ 20 km s -1 in a dBarth et al.1 
l2002t iCappella ri fc Emsellem| |2f )f)4D . F or the SAURON 
galaxy sample. lEmsellcm et al.l (120041 ) found scatter of 
15-18 kms -1 in comparisons between their pPXF-based 
measurements of a and previous a values in the litera¬ 
ture. For our trials, we adopt the conservative assump¬ 
tion that v and a each have errors of 20 km s -1 . We have 
run pPXF using an alternative subset of MILES stellar 
templates and found that the resulting variations in v 


and a were well within this range. 

In order to assess the kinematic error terms for each 
line index, we have tested five spectra from NGC 1023. 
The innermost spectrum spans radii from 0.5" to 2.0", 
with S/N = 119 over the wavelength range used for the 
pPXF fit, and a = 216kms -1 . The outermost spectrum 
spans 50"-75" and has S/N = 38 and er = 94kms -1 . For 
a given spectrum, we perform two sets of trials wherein 
we vary the assumed velocity and velocity dispersion, vo 
and (Jq. In one set of trials, we fix a o to the best-fit 
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a from pPXF and sample vq from a normal distribution 
with a dispersion of 20kms~ 1 (centered on the best-fit v 
from pPXF). For each of 100 trials sampling vq , we con¬ 
volve the spectrum from <to to a = 230 km s, shift to 
Arest = A 0 bs/ (1 + vq / c ), and measure the absorption line 
indices. For each absorption feature, we define a system¬ 
atic error term e v o as the standard deviation in index 
strength over all 100 trials. Our second set of trials as¬ 
sesses an analogous error term e a Q by employing a nearly 
identical procedure whereby v$ is held fixed and cto is 
sampled from a normal distribution with a dispersion of 


20kms _1 . 

In Table [3] we present the values of e„o and e CT o for 13 
different line indices, averaged over our five test spec¬ 
tra. We make particular note of the Nal index, which 
is 2-3 times stronger at the galaxy center than at large 
radii. This index exhibits a similar absolute value of 
e„o ~ 0.025 A at all radii, while e a o scales approximately 
linearly with the index strength. The trends in e„o and 
e CT o for both galaxies’ Nal indices are illustrated in Fig¬ 
ure [TGl Figures 1X71 and [TS1 show individual spectra of the 
Nal region for NGC 1023, with variations in vq and a 0 
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Figure 17. (a) line strength of Nal feature vs. input systemic 

velocity (no), for a spectrum covering major-axis radii from +0.5" 
to +2.0" in NGC 1023. (b) Nal line strength vs. input velocity 
dispersion (cro) - In each case, the spectrum is convolved from <ro 
to cr = 230 km s^ 1 before measuring the line index, (c) rest-frame 
spectrum of the Nal feature, before convolution to o = 230 kms -1 . 
The red segment indicates the Nal line region, and the blue seg¬ 
ments indicate the pseudo-continuum regions. The blue dashed 
line is the linear fit to the pseudo-continuum, used in the index 
measurement, (d) rest-frame spectrum at two separate velocity 
shifts, (e) rest-frame spectrum convolved to 230 kms -1 from two 
separate op values. 


Figure 18. Panels are the same as Figure [TT] This spectrum 
spans major-axis radii from —50" to —30" in NGC 1023. The ve¬ 
locity dispersion is smaller than the central spectrum in Figure [TT] 
and the Nal line is significantly weaker. 
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and the resulting smoothed spectra. Although the inner 
spectrum (Figure fTTl) has a much deeper Nal feature and 
much larger intrinsic er, the change in index strength with 
respect to vq is similar to the outer spectrum (Figure flSl) . 

As indicated in Table [3j systematic errors from kine¬ 
matic fitting are typically ~ l%-5% of the line index 
strength. In cases of low statistical noise, the e„o and/or 
Co-o term is comparable to the level of random variance 
between different data subsets. This is true for index 
measurements based on high-S/N spectra near the cen¬ 
ter of each galaxy, and even toward large radii for the Fe, 
Mg b, and Call features. On the other hand, the Nal 
and MgI0.88 features are relatively sensitive to the ve¬ 
locity shift and broadening of the underlying spectrum, 
such that the kinematic error terms make a substantial 
contribution to the overall error budget even when S/N 
is modest. 

A.2: Contamination from Emission Lines 

Emission lines from warm gas are visible contami¬ 
nants in our spectra throughout NGC 2974 and at the 
center of NGC 1023. We have experimented with two 
well-known routines to remove the emission component 
and ass ess the under lying stellar abs orption features: 
pPXF bv lCappellari fc Emsellcml (I2004T) and GANDALF by 
iSarzi et al.l (1200(f ). In addition to the stellar template li¬ 
brary used to fit a galaxy spectrum, pPXF supports Gaus¬ 
sian emission line templates at the rest wavelengths of 
several gas species. The strength of each emission line is 
varied freely, while the emission line kinematics (v and 
a) are assumed to match the stellar kinematics. Using 
the output template data provided by pPXF, we separate 
the stellar template and emission line components, and 
subtract the sum of the emission lines from our origi¬ 
nal galaxy spectrum. This is performed in the same run 
where we measure v and a. 

The GANDALF routine adds more flexibility to the 
emission-line fitting: the gas component(s) are permit¬ 
ted to have different kinematics from the stars, and the 
relative flux and kinematics of each emission line may ei¬ 
ther be fit freely or coupled to another line. In practice, 
the routine performs best when the stellar kinematics are 
first estimated using pPXF (while masking the emission 
lines or fixing their kinematics as above), and supplied as 
an initial guess for the subsequent GANDALF fit. We have 
used GANDALF successfully for NGC 2974, whose strong 
emission lines are easily recognized. As before, we isolate 
the emission-line component of the best-fitting spectrum 
and subtract it from our original galaxy spectrum. Any 
emission in NGC 1023 is too subtle for GANDALF to sepa¬ 
rate cleanly, and it confuses the low-order stellar contin¬ 
uum with a superposition of extremely broad emission 
lines. Therefore, we restrict our use of GANDALF to NGC 
2974. 

We find that the best-fitting emission spectrum varies 
for different sets of stellar templates used in pPXF, and 
with the degree to which the relative emission strengths 
of different species are allowed to vary in GANDALF. We 
therefore define a systematic error term e gas as the stan¬ 
dard deviation of trial measurements for a given absorp¬ 
tion line index, when different settings are used for the 
emission line removal. For NGC 2974, the trials include 
two alternative stellar template libraries - each a subset 
of the empirical MILES library - for pPXF, and two set- 
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Figure 19. Galaxy spectra before and after emission line correc¬ 
tion, zoomed in on the H (3 and [OIII] features. In each case the 
thin grey line represents the uncorrected spectrum, and the thick 
colored line represents the spectrum after correction with emission 
line templates in pPXF. From top to bottom: central bin of NGC 
1023; central bin of NGC 2974; outermost bin for NGC 2974. 

tings for the relative emission line strengths in GANDALF. 
In one GANDALF trial the flux of each emission line is 
treated as a free parameter, and in the other trial we 
fix the relative strengths within the [Oil], [OIII], [NI], 
and [ArlV] multiplets. In NGC 2974, Balmer emission 
lines impact the CNi, bTiO, and H/3 absorption features. 
For H/3 we find that e gas ss 6% of the absorption line 
strength, and bTiO is impacted at the « 5% level as a 
result of H/3 contaminating the red pseudo-continuum. 
The Mg b and C24668 indices exhibit small systematic 
errors from [NI] and [ArlV] contamination, respectively. 
For the average index measurements presented in ([5] and 
Figures OH] we have adopted pPXF as our default tool for 
emission line removal, as we find it performs better at 
distinguishing between relatively narrow emission lines 
and broad variations in the galaxy continuum relative to 
the stellar templates. 

For NGC 1023 we determine e gas using two trials: a 
pPXF trial with our default template library, and a trial 
with no emission line fitting. A third trial with an alter¬ 
native template library for pPXF does not differ substan¬ 
tially for NGC 1023, and we exclude it so as not to dilute 
e gas . As with NGC 2974, we find a ~ 6% impact on the 
H/3 index, and possible low levels of contamination for 
the CNi and Mg b absorption features. In Figure [12] 
we show example spectra from both galaxies, before and 
after emission line removal. 

A.3: Overlap of redshifted Nal and telluric H 2 O 

We have attempted to remove telluric absorption fea¬ 
tures from our galaxy spectra by comparing them to 
a grid of model atmospheric transmission spectra over 
the observed wavelength range 9310-9370 A. Although 
our corrected spectra look reasonable by eye, we aim to 
quantify the level of uncertainty near the 8190A Nal 
feature and the corresponding bias in our measurement 
of the Nal line index. To this end, we have performed 
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Figure 20. Galaxy spectra near the Nal absorption feature, after varying the strength of telluric absorption bands. Left: spectra from 
our central and outermost bins in NGC 1023. Right: spectra from our central and outermost bins in NGC 2974. 
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Figure 21. Comparison of Nal measurements in NGC 1023, for various levels of telluric correction. The rotation curve of the galaxy 
will shift stellar absorption features with respect to imperfectly corrected telluric absorption bands, driving asymmetries in the line index 
measurements from opposite sides of the galaxy. Left: Nal measurements after weakening our best-fitting atmospheric transmission spectra 
by 10%. Middle: Nal measurements using our best-fitting transmission spectra. Right: Nal measurements after strengthening our best¬ 
fitting transmission spectra by 10%. The middle and right panels both represent a plausible strength for telluric absorption features in our 
spectra, whereas the left panel exhibits clear asymmetry between opposite sides of the galaxy as the result of poor telluric correction. 


trials where we multiplied the strength of our best-fitting 
transmission spectrum for each science exposure by val¬ 
ues of 0.9, 1.1, and 1.2, and repeated all data processing 
and analysis steps starting from telluric division of in¬ 
dividual frames. In Figure [20] we display the resulting 
spectra for each galaxy, from our central spatial bin and 
an outer spatial bin. 

For NGC 1023 (Figure l20k ). contamination by the deep 
H 2 O feature at 8230 A is visible for telluric correction 
of -10% and +20% relative to the best fit. Beyond a 
crude visual estimate of our possible error range for tel¬ 
luric correction, we can take advantage of the fact that 
NGC 1023 and NGC 2974 both have rotational veloci¬ 
ties ~ 200kms , such that the Nal feature is shifted 


by ~ 10 A over the length of our slit. If the Nal line or 
pseudo-continuum regions are severely contaminated by 
a telluric feature, we should see systematic differences 
in our index measurements from opposite sides of the 
galaxy. This is illustrated in Figure I2T1 For each panel, 
we have overplotted all four subsamples used to compute 
our random errors in the Nal index (see H4.1D . For the 
trial where we adjusted our transmission spectra by -10% 
(Figure HB l ), our measurements of the Nal index on op¬ 
posite sides of the galaxy are clearly offset. We see similar 
asymmetry for the trial with +20% adjustment, though 
it is not displayed in Figure [2TJ Our other two trials (no 
adjustment and +10%, corresponding to Figures 12Tb and 
EDO both exhibit overlapping measurements from oppo- 
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site sides of the galaxy and do not produce obviously poor 
spectra. We therefore estimate our plausible systematic 
error from telluric correction in terms of the difference 
between our Nal index measurements from these latter 
two trials. Specifically, we have adjusted our final Nal 
index measurement in each spatial bin of NGC 1023 to 
the average value from the two trials. In each bin we 
define the 1 -a systematic error term e te i as the percent 
deviation between the two trials, equal to -4^ times the 
percent difference in index measurements. 

Our random errors in each line index are derived from 
subsets of data spanning both sides of the galaxy f M.ll) . 
and therefore some of the scatter in Figure I2T1 is already 
incorporated in the random error term for Nal. How¬ 
ever, Figure [21] also exhibits a systematic offset between 
the three panels, indicating that poor telluric correction 
introduces an overall bias in addition to increased scat¬ 
ter from the shift in observed wavelengths of Nal. We 
therefore keep the etei term as defined above, to account 
for the plausible bias level from over- or under-corrected 
telluric features. 

We observed NGC 2974 during the second half of our 
observing night, and telluric H 2 O in our spectra is not as 
strong as for NGC 1023 (see Figure 120b j. Yet our -10% 
and +20% atmospheric transmission trials once again 
produce offsets in the Nal index measurements from op¬ 
posite sides of the galaxy. Therefore we define e te 1 for 



Figure 22. Comparison of the Nal and Nalgugg indices defined 
in Table[2] In each panel the red segment indicates the line region, 
and the blue segments indicate the pseudo-continuum regions. The 
blue dashed line is the linear fit to the pseudo-continuum, used in 
the index m easu rement. The example galaxy spectrum is the same 
as in Figure [771 prior to convolution to <r = 230kms _1 . 


NGC 2974 in the same manner as for NGC 1023, and 
adjust our average index values accordingly. For NGC 
2974 we compute a maximum e te i of 3.8%, versus 13.4% 
for NGC 1023. For both galaxies we add et e i in quadra¬ 
ture with other systematic terms, as indicated in Tablc[3] 

B: ALTERNATIVE DEFINITIONS OF THE SODIUM 8190A 
AND FEH WING-FORD INDICES 

In rec ent inve stig ations of IMF-sen s itive absorption 
features, lCvD12l and lLa Barbera et all (|2013f l have pro¬ 
posed two different definitions of the Nal absorption in¬ 
dex near 8190 A. We compare the corresponding line 
and pse udo-con tinuum regions in Figure 1221 In the ver¬ 
sion by [C vD12t both pseudo-continuum regions are im¬ 
mediately adjacent to the line region, such that the red 
pseudo-continuum lies directly on the overlapping TiO 
band. The NalsDSS index from lLa Barbera et all ( 2013 ) 
defines a narrower line region and moves both pseudo- 


0.0 0.5 1.0 f (k f.5 ) 2.0 2.5 




V I X eff 

Figure 23. Ratio of Nal and NalsDSS to (Fe) as a function 
of radius. Left : NGC 1023, with the outermost spatial bin ex¬ 
tending to 1.6r e ff. Right: NGC 2974, with the outermost spa¬ 
tial bin extending to 0.8r e ff. The grey shaded area in each 
panel indicates the sca l e of t he single-aperture measurement by 
Ivan Dokkum fc Conrovl (120121) . 
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Figure 24. Comparison of the FeH index definition from ICvD12l and three variants explored herein. In each panel the red segment 
indicates the line region, and the blue segments indicate the pseudo-continuum regions. The blue dashed line is the linear fit to the 
pseudo-continuum, used in the index measurement. The example galaxy spectrum is the innermost bin (—5" to +5 // ) for NGC 1023, 
with no velocity convolution. For visualization purposes, the spectrum has been divided by a fourth-order polynomial approximating the 
continuum from 9800-10,050 A. 
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Figure 25. Ratio of FeH to (Fe) as a function of radius, using the four definitions of the FeH index in Figure [24l Left: NGC 1023, with 
good data for FeH extending to 1.0r e ff. Right: NGC 2974, wit h g ood data fo r F eH extendin g to 0.26r e fj. The grey shaded area in each 
panel indicates the scale of the single-aperture measurement bv lvan Dokkum K: Conrovl 12012 1*1. 


continua to regions outside the Nal+TiO blend. This 
generates a shallower linear fit to the continuum and pro¬ 
duces larger equivalent widths. However, we find that 
the offset between our measured Nal and NalsDSS line 


strengths is nearly the same at all radii in NGC 1023 and 
NGC 2974. In Figurel251we show the ratios Nal/(Fe) and 
NaIsDSs/(Fe) as a function of radius in NGC 1023 and 
NGC 2974. For both definitions, the relative strength of 
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the sodium index rises steeply from 0.1 r e Q to the galaxy 
center. For all figures and di scussion above, we follow 
the Nal index definition from lCvD12l 
The relative strengths of the FeH, Nal, and (Fe) ab¬ 
sorption features are crucial for our interpretation of 
abundance versus IMF gradients in NGC 1023 and NGC 
2974. The FeH index is notoriously difficult to mea¬ 
sure, as spectra near 9900 A suffer from bright telluric 
emission lines and relatively low instrumental through¬ 
put. While the Wing-Ford feature corresponds to the 
principal FeH bandhead in this region, the shape of the 
nearby pseudo-continuum is also influenc ed by overlap¬ 
ping bands of Til, TiO, and CrH (e.g.. IMcLean et all 
120031 : [Cushing etahl I2005L CvD12). To check whether 
our observed radial trends in FeH are sensitive to the pre¬ 
cise placement of the line and pseudo-continuum regions, 
we have repeated our measurements using four variants 
of the index definition. These are illu strated in Figure [Ml 
The first variant is the definition bv ICvDlH which we list 
in Tableland use for all analysis and figures above. In 
this version, the line region is defined to match the deep¬ 
est part of the FeH bandhead, possibly reducing contami¬ 
nation from overlapping TiO. Our second variant instead 
uses a larger portion of the FeH trough (9902.3-9962.3 
A) and moves the red pseudo-continuum accordingly to 
9967.3-9982.3 A. The third and fourth variants match the 
first two, except with the blue pseudo-continuum moved 
from 9852.3-9877.3 A to 9862.3-9887.3 A. 

We display the radial trends of FeH/(Fe) for each vari¬ 
ant in Figure [25] None of the adopted FeH variants 
alter our main finding: that FeH decreases relative to 
(Fe) going from 0.1r e ff into the center of each galaxy, 
and thereby opposes the trend in NaI/(Fe) at these 
same radii. Yet there is evidence of some deviation at 
r > 0.5r e ff in NGC 1023, such that the variants using 
the full absorption trough (variants 2 and 4) show an 
increase i n FeH to ward large r, whereas the original def¬ 
inition by ICv Dll and our similar variant 3 show a slight 
decrease toward large r. This may be caused by a faint 
sky line near 9970 A rest, which overlaps with the red 
pseudo-continuum in our full-trough variants. Traces of 
this line are visible in Figure [3}) and for the two outer¬ 
most spectra in Figure 0]:. 
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